Exotic phases induced by strong spin-orbit coupling in ordered double perovskites 
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We construct and analyze a microscopic model for insulating rock salt ordered double perovskites, with the 
chemical formula A2BB'06, where the B' atom has a Ad} or 5d^ electronic configuration and forms a face 
centered cubic (fee) lattice. The combination of the triply-degenerate t2g orbital and strong spin-orbit coupling 
forms local quadruplets with an effective spin moment j = 3/2. Moreover, due to strongly orbital-dependent 
exchange, the effective spins have substantial biquadratic and bicubic interactions (fourth and sixth order in 
the spins, respectively). This leads, at the mean field level, to three main phases: an unusual antiferromag- 
net with dominant octupolar order, a ferromagnetic phase with magnetization along the [110] direction, and a 
non-magnetic but quadrupolar ordered phase, which is stabilized by thermal fluctuations and intermediate tem- 
peratures. All these phases have a two sublattice structure described by the ordering wavevector Q = 27r(001). 
We consider quantum fluctuations and argue that in the regime of dominant antiferromagnetic exchange, a non- 
magnetic valence bond solid or quantum spin liquid state may be favored instead. Candidate quantum spin 
liquid states and their basic properties are described. We also address the effect of single-site anisotropy driven 
by lattice distortions. Existing and possible future experiments are discussed in light of these results. 



PACS numbers: 71.70.Ej,71.70.Gm,75.10.-b 

I. INTRODUCTION 

In magnetic Mott insulators with quenched orbital degrees 
of freedom, weak spin-orbit coupling (SOC) only leads to a 
small correction to the usual spin exchange Hamiltonian in 
the form of single-site anisotropy and Dzyaloshinskii-Moriya 
interactions.^^ In the presence of strong SOC, however, a 
completely different physical picture emerges, in which spin 
itself is not a good quantum number, and magnetic anisotropy 
is usually large. Generally, strong SOC is common in the Lan- 
thanides, in which the relevant 4f-electrons are very tightly 
bound to the nucleus. The tight binding shields the electrons 
from crystal fields, which tend to split the orbital degeneracies 
involved in SOC, and moreover reduces exchange, which also 
competes with SOC. 

While more rare, strong SOC is becoming an increasing 

focus in d-electron systems, in which electrons are more delo- 

calized than in the Lanthanides, and more diverse phenomena 

can be expected. For instance, strong SOC can be expected 

in 5d transition metal compounds, which have large intrinsic 

atomic SOC due to their high atomic weight. In this category, 

many Ir-based magnets have been studied recently^ Lighter 

transition metals may also exhibit strong SOC if competing 

effects such as crystal fields and exchange are suppressed, 

e.g. by choosing crystal structures with high-symmetry and 

well- separated magnetic ions, respectively. An example of 

this type is the "spin-orbital liquid" state observed in the Fe- 

based spinel FeSc2S4,^'^ which is believed to be driven by 
3OC.^M0j 



In this paper, we consider the case of insulating magnetic 
ordered double perovskites. Structurally, ordered double per- 
ovskites (with the chemical formula A2BB'06) are derived 
from the usual perovskites ABO3 by selectively replacing half 
the B ions with another species, denoted B'. We focus on the 
case in which the B ions are non-magnetic and the B ' ones are 
magnetic. Because of the difference in the valence charges 
and ionic radius between B and B' ions, the magnetic B' ions 
form an fee lattice structure with a lattice constant double of 
the original cubic one. Many ordered double perovskites in- 
corporate strong intrinsic SOC, as B' ions are commonly 4d 
and 5d transition metals. Moreover, the large B'-B' separation 
weakens exchange, similarly to FeSc2S4. Here, we construct 
an appropriate microscopic model for the most quantum of 
these materials (a list may be found in Table |ll]), in which the 
magnetic ion contains a single unpaired S = 1/2 spin. 

The physics is strongly influenced by the combination of 
the orbital degeneracy of the t2g multiplet, which acts as an ef- 
fective ^ = 1 orbital angular moment. Due to strong SOC, this 
combines with the 5* = 1/2 spin to induce an effective total 
angular momentum j = 3/2 description of the system. More- 
over, due to the orbitally-dependent exchange, the interaction 
of these j = 3/2 contains large biquadratic (fourth order in 
spin operators) and triquadratic (sixth order in spin operators) 
interactions. These support exotic phases not easily found in 
systems with dominant bilinear spin exchange. 

Analysis of the microscopic model shows that the strong 
SOC enhances quantum fluctuations and leads to several in- 
teresting phases: (1) an unconventional antiferromagnet (de- 
noted AFM) in which the magnetic octupole and quadrupole 
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moments rather than the dipole moment are dominant, (2) 
an unusual non-colHnear ferromagnet (denoted FMllO) with 
a doubled unit cell and magnetization along the [110] axis, 

(3) a (biaxial) "spin nematic" phase with quadrupolar order 
but unbroken time reversal symmetry and, more speculatively, 

(4) a possible quantum spin liquid (QSL) phase. Phases (1), 
(2), and (4) are low temperature phases and persist as ground 
states, while the spin nematic, phase (3), occurs in a broad 
intermediate temperature range below the paramagnetic state 
but above any magnetic ordering temperature. 
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FIG. 1. (Color online) Mean field T = phase diagram for the 
model Hamiltonian in Eq. ( |27| l. AFM denotes the "antiferromag- 
netic" ground state given by Eq. ( [39] l and Eq. ( |40] ), FMllO de- 
notes the ferromagnetic ground state with an easy axis oriented along 
[110], given by Eq. |5T}. The FMIOO state, which is ferromagnetic 
with easy axis along [100] appears in the narrow band between the 
two phase phase boundaries. In the figure, J = 1. 

States with magnetic multipole order are more often ob- 
served in f-electron systems where crystal field effects are less 
important than SOC.^^ As a consequence, the atomic wave- 
functions are total angular momentum eigenstates, in which 
the spin and orbital degrees of freedom are highly entangled. 
This leads to highly non-Heisenberg exchange between the lo- 
cal moments, which is described by interaction of higher mag- 
netic multipole operators. Such interactions may drive multi- 
polar order, as suggested for instance in URu2Si2.^^ Recently 
this has been suggested to also occur in d electron systems 
with unquenched orbital degeneracy and sufficient SOC.^^ We 
find a similar mechanism at work in the AFM phase. 

A ferromagnetic state is not in itself unusual, though such 
is relatively uncommon in insulators. However, cubic ferro- 
magnets with an easy axis oriented along the [110] direction 
is quite uncommon. This can be understood from the Lan- 
dau theory for a ferromagnet: the usual fourth order cubic 
anisotropy term favors either [100] or [111] orientation, de- 
pending upon its sign, but never [110]. To obtain a [110] easy 
axis, one requires sixth order or higher terms to be substantial, 
making this rare indeed. Remarkably, such [110] anisotropy 
has been observed in experiments on Ba2NaOs06.^'^ 

Both the above states, when heated above their magnetic 
ordering temperatures, allow on symmetry grounds for an in- 
termediate phase which is time-reversal symmetric but with 



quadrupolar order - the spin nematic. Applying the mean field 
theory at T > 0, we indeed find such a phase in a broad range 
of parameter space. While spin nematic states have been sug- 
gested previously in NiGa2S4^^ it has not been established 
in that material. The mechanism for quadrupolar order here is 
much more transparent and robust than in that case. 

The above three phases, while somewhat unconventional, 
may be obtained within a mean-field analysis. A QSL state, 
however, cannot be described by any mean field theory, and 
is considerably more exotic. The search for a QSL, which is 
a state in which quantum fluctuations prevent spins from or- 
dering even at zero temperature, is a long-standing problem 
in fundamental physics. Since the possibility of a QSL was 
suggested by Anderson in the early 1970s,'^this has been an 
active area for theory and experiment. Despite the current ma- 
turity of the theory for QSL,^^ the experimental confirmation 
of the existence of such an exotic phase is still elusive. Very 
commonly geometrical frustration is thought to be a driving 
mechanism for QSL formation, and consequently most re- 
search (both theoretically and experimentally) has been de- 
voted to systems of this type, such as triangular,^^ kagome,«^ 
hyperkagome^ and pyrochlore lattices. 

Here we suggest a different route, in which quantum fluctu- 
ations are enhanced primarily by strong SOC, rather than geo- 
metrical frustration. In fact, the magnetic ions in ordered dou- 
ble perovskites reside on a face centered cubic (fee) sublattice, 
which can be viewed as edge- sharing tetrahedra, and is some- 
what geometrically frustrated. Without strong SOC, however, 
this frustration is weak, and indeed the classical Heisenberg 
antiferromagnet on the fee lattice is known to magnetically or- 
der into a state with the ordering wavevector 27r(001).'2SlThe 
tendency of the simple fee antiferromagnet to order may be 
partially attributed to its large coordination number (z = 12), 
which leads to mean-field like behavior. By contrast, strong 
SOC induces effective exchange interactions very different 
from Heisenberg type, with strong directional dependence that 
may make a QSL more favorable. To make this suggestion 
more concrete, we propose a natural wavefunction for a QSL 
in our model, and discuss the physical properties of such a 
state. 

We now outline the main results of the paper, and how they 
are presented in the following sections. In Sec. |ll| we show 
that strong SOC leads to an effective j = 3/2 local moment 
on each B' site. We write down a model Hamiltonian which 
includes three interactions: nearest neighbor (NN) antiferro- 
magnetic (AFM) exchange, J, NN ferromagnetic (FM) ex- 
change, J^ and electric quadrupolar interaction, V. These in- 
teractions are all projected down to the effective j = 3/2 man- 
ifold, which induces many terms beyond the usual quadratic 
exchange. Indeed, because of the four-dimensional basis of 
spin-3/2 states, the resulting Hamiltonian can be thought of as 
an anisotropic F matrix model. We then discuss the symme- 
try properties of the projected Hamiltonian. Surprisingly, we 
find that, in the limit of vanishing FM exchange, the Hamilto- 
nian has a "hidden" global SU(2) symmetry despite its com- 
plicated appearance. 

In Sec.|Ill| we consider the mean field ground states of the 
model, characterized by local (single-site) order parameters. 
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In Sec. Ill A we begin by considering the more accessible limit 
in which strong uniaxial single- site anisotropy (due e.g. to a 
tetragonal distortion of the crystal) lifts the "orbital" four-fold 
degeneracy of j = 3/2 quadruplets down to easy-axis or easy- 
plane Kramer's doublets. In these limits, the effective Hamil- 
tonian in the reduced phase space is mapped onto that of an 
XXZ antiferromagnet which can be understood even without 
mean field theory. Next, in Sec. IIIB we carry out T = 
mean field theory for the case of cubic symmetry. Here we 
find the AFM state and two ferromagnetic states (the FMllO 
state and another state with a [100] easy axis). The T = 
mean field phase diagram is shown in Fig. [T] Finally, having 
described the situations with strong and vanishing single- site 



anisotropy, we determine in Sec. |III C| the mean-field phase 
diagram for intermediate anisotropy. 



In Sec. IV we identify the multipolar order parameters of 
the three ordered phases and analyze the T > behavior by 
mean-field theory. Here we find the quadrupolar phase, and 
discuss several phase transitions which occur. We also dis- 
cuss the behavior of the magnetic susceptibility in different 
parameter regimes. 

In Sec. |Vj we consider quantum effects beyond the mean 
field theory. First, we carry out a spin wave calculation, which 
determines the collective mode structure, and also shows that 
in the regime where nearest-neighbor antiferromagnetic ex- 
change is dominant (small J' and V) quantum fluctuations are 
large and may destabilize the ordered AFM phase. Therefore, 
we consider possible non-magnetic ground states, both of va- 
lence bond solid and quantum spin liquid (QSL) type. We 
formulate a slave-fermion theory with four-component spin 
s = 3/2 fermions, in such a way that at mean field level the 
hidden SU(2) symmetry is preserved and the correct ground 
state, the analog of a singlet in the usual Heisenberg model, is 
obtained for a single pair of nearest-neighbor sites. The cor- 
responding mean field theory naturally includes the intrinsic 
spatial anisotropy of the strong SOC limit. We analyze two 
different mean field ansatze, with zero and 7r-flux. In both 
cases the mean field Hamiltonian respects all the symmetries 
of the original spin Hamiltonian. The 7r-flux state is found to 
have lower mean field energy. For both states, the spinons are 
at quarter-filling, leading to a spinon Fermi sea. There is no 
Fermi surface nesting and we expect that this spinon Fermi 
surface should be stable against weak perturbations. Predic- 
tions based on the picture of spinon Fermi surface are made. 



Finally in Sec. |Vlj we compare our theoretical prediction 
with current experimental findings and suggest further direc- 
tions for theory and experiment. 



II. MODEL AND SYMMETRY 

A. Spin-orbit interaction and hybridization of atomic orbitals 

The magnetic ions B' (Os^+, Re^+, Mo^+) found in the 
ordered double perovskites in Table [ll| all have one electron 
in the triply degenerate t2g multiplet. The atomic spin-orbit 



interaction projected down to the t2g triplet is written as 

Hso = -A / • ^ , (1) 

in which the total angular momentum quantum numbers of 
these operators are / = 1,6' = 1/2. The effective orbital 
angular momentum / comes from the projection of orbital an- 
gular momentum L onto the t2g triplets. 



-I 



(2) 



Here Vt^, = T.a=yz,xz,xy is the projection operator 

to the t2g manifold. The eigenstates of P with eigenvalues 
m = 0, ±1 and with eigenvalues ct = ±1/2 =t, J,, written 
in terms of the usual t2g states are 



|0,ct) = |<J; |±1,ct) = 



(3) 



This interaction favors j = 3/2 (j = I + S) quadruplets 
over j = 1/2 doublets by an energy separation 3 A/2. In the 
strong spin-orbit interaction limit, the local Hilbert space is 
restricted to four low-lying states 



\da) = ^C^^|m,(j), 



where a = ±3/2, ±1/2 is the eigenvalue and 



(4) 



-, a 



(5) 

is a Clebsch-Gordan coefficient. In the materials under con- 
sideration, A is indeed a very large energy scale (fraction of 
an eV), justifying the strong SOC limit. 

Every operator expressed in terms of spin and orbitals must 
therefore be projected into this subspace and its projection can 
be written in terms of j = 3/2 angular momentum operator. 
For example. 



Vs I Vs 



(6) 
(7) 



Here Vs is the projection operator into the j = 3/2 quadru- 
plets. Furthermore, for the magnetic moment M for electrons 
in atomic d orbitals, we have 



M = V3[2S^{-l)]V3 = 0. 



(8) 



The vanishing magnetic moment is quite remarkable and par- 
tially explains why the compounds have small magnetic mo- 
ments in comparison with spin- ^ systems without orbital de- 
generacy. 

In reality, the measured magnetic moments are nonzero be- 
cause the atomic 4d or 5d orbitals strongly hybridize with p or- 
bitals at the oxygen sites that form an octahedron surrounding 
each B' site. For instance, for Ba2NaOs06, the hybridization 
energy is estimated to be of the order of electron volts and 
comparable to the energy gap between Os d and O p states. 
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For this reason, it is more appropriate to think in terms of 
molecular orbitals with mixed d and p character. For exam- 
ple, molecular xy orbitals are written as 



W ) 



(9) 



where \d^^^y) is the state corresponding to one electron in the 
xy orbital and spin cr on site z, and \p^^j.y) is a linear combi- 
nation (with xy symmetry) of states that have a singlet on the 
dxy orbital and one hole on an oxygen site 



) + |Pr-e.,J + |P^ 



(10) 

where e^^y are real space vectors from the B' site to neighbor- 
ing oxygens along x or y directions. The mixing parameter r 
is of order tdp/ A, where tdp is the hopping matrix element be- 
tween d and p orbitals and A is the gap to oxygen p states. In 
the limit of strong spin-orbit interaction, we must project into 
four low-lying molecular orbitals which are a superposition of 
the four atomic states with j = 3/2 and p states 



I A, a) — ^ \Di,m) • 



(11) 



where the sum is over nearest neighbor sites in the XY planes, 
and the corresponding terms for YZ and XZ planes are ob- 
tained by the obvious cubic permutation. Here the operators 
Si^xy and rii^xy denote the spin residing on xy orbital and or- 
bital occupation number at site i, respectively. In terms of spin 
and orbital angular momentum operators acting on site i. 



(14) 
(15) 



Throughout this paper, we use the subindices (z, xy) to de- 
note the site and orbitals, superindex (/i = x^y^z) to denote 
the spin component, and capital letters (XY, XZ, YZ) to de- 
note the planes. With these definitions, we note that the single 
occupancy condition at each site, which defines the Mott in- 
sulating state, becomes 



^xy H~ ^i,xz H~ ^i,yz — !• 



(16) 



Moreover, from Eq. ([T5j, orbitally-resolved spins satisfy 



^i,xy H~ ^i,yz H~ ^i,xz — 



(17) 



While the atomic magnetic moment in Eq. ([8| vanishes, there 
is a nonzero contribution to the molecular M from holes in p 
orbitals. After taking the projection into j = 3/2 states, the 
coupling of the molecular orbital to a magnetic field reads 



Hz = -giiBh' j, 



(12) 



where g = r^/[3(l + r^)] is the Lande factor, and /i^ is the 
Bohr magneton. 




B. Exchange interactions and electric quadrupolar interaction 

In the last subsection, we discussed the effect of strong 
spin-orbit interaction in determining the local degrees of free- 
dom and pointed out that every operator must be projected 
into the j = 3/2 quadruplets. In this subsection, we introduce 
the interactions between the local moments, and discuss the 
mechanics of the projection. 

The first interaction to consider is nearest-neighbor anti- 
ferromagnetic exchange, through the virtual transfer of elec- 
trons through intermediate oxygen p orbitals. These processes 
are strongly restricted by symmetry. For example, in XY 
planes, only electrons residing on d^y orbitals can virtually 
hop to neighboring sites via Px and Py orbitals of the interme- 
diate oxygen sites. The exchange path and relevant orbitals 
are depicted in Fig. [2j Alternatively, one can interpret this 
process as kinetic exchange between molecular Dxy orbitals, 
which are mixtures of the transition metal d state and p states 
on the neighboring four oxygens (see Eq. (|9])). As a con- 
sequence, the antiferromagnetic exchange interaction can be 
written Hex-i = HfJi + n^^_^ + Hf^^i, where 



XY 



1 = ^ E 

{ij)exY 



^i,xy ' 



1 



(13) 



FIG. 2. (Color online) Left graph: The NN AFM exchange path (B'- 
O-O-B^; right graph: The NN FM exchange path with intermediate 
orthogonal p orbitals at O sites. 

The second interaction is the nearest-neighbor ferromag- 
netic exchange interaction. This interaction is due to the spin 
transfer through orthogonal orbitals at the intermediate oxy- 
gen sites in the exchange path, as shown in Fig. |2] For two 
sites j in the XY plane, this ferromagnetic exchange is writ- 
ten as 



K 



XY 



ex-2,ij ^ ~J [^i,xy ' {^j,yz + ^j,xz) + ^ j)] 

3J' 

-—r- [ni^xyirij^xz + rij^yz) + ^ j)] , (18) 



where the xy orbital only interacts with yz and xz orbitals at 
neighboring sites. Applying the single-occupancy constraint, 
the nearest-neighbor ferromagnetic exchange interaction can 
be simplified, up to a constant, to 

^&^-2 = — J' ^ [^i,xy ' {^j,yz + ^j,xz) + ^ j)] 

{ij)exY 
3J' ^ 



2 / ^ ^i,xy^j,xy • 



(19) 
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Microscopically, J' / J ^ 0{Jh /Up) where Jh and Up are 
the Hund's coupling and Hubbard Coulomb interaction at the 
oxygen site, respectively. 

The third interaction is the electric quadrupole-quadrupole 
interaction. The 4d or 5d electron carries an electric 
quadrupole moment, and the interaction between these mo- 
ments may not be negligible because of the long spatial ex- 
tent of the molecular orbitals. Calculating the direct electro- 
static energy between all possible orbital configurations for 
two electrons residing in neighboring sites in an XY plane, 
we obtain the quadrupole-quadrupole interaction 

1-1^^ — Vn ■ n 



V 

2" 

13V 



(y^i,yz^j,yz ~l~ ^i,xz^j,xz) 



12 
19V 

H {'^i,yz'^j,xz ~^ '^i,xz'^j,yz) ' (20) 

Here V > is defined as the Coulomb repulsion between two 
nearest-neighbor xy orbitals on XY planes. If Q is the magni- 
tude of the electric quadrupole and a is the lattice constant of 
the fee lattice, we have V = 9\/2Q^ j . In general, the main 
contribution to Q comes from the charge at the oxygen sites, 
hence the larger the hybridization, the larger the value of V . 
Using the single-occupancy constraint in Eq. ([16]) and sum- 
ming over sites, the quadrupole-quadrupole interaction sim- 
plifies to 



O/XY 
'^quad 



E 

(ij)GXY 

W 



0K> 



(21) 



in which we have ignored an unimportant constant. 

The minimal Hamiltonian for the cubic system contains all 
three of these exchange interactions in addition to the on-site 

soc. 



^ — ^ex-1 + + ^( 



quad 



(22) 



Since we are interested in the limit of strong spin-orbit in- 
teraction, we need to project H onto the j = 3/2 quadruplets 
at every site. As an example, we write down the projection for 



^i,xy and Tii^xy-) 



nx 

i,xy 




AZ ^X 'Z 

2 Ji Ji Ji 


(23) 


oy _ 

i,xy 


1 -y 


-z -y -z 

2 Ji Ji Ji 


(24) 


^i^xy 




-z -z -z 

2 Ji Ji Ji 


(25) 


'^i,xy 


3 

4 ~ 




(26) 



in which, O = V3 O Vs. Spin and occupation number op- 
erators for other orbitals can be readily generated by a cubic 



permutation. After the projection, the minimal Hamiltonian 
reduces, up to a constant, to 



^ — ^ex-l + + ^( 



quad 



(27) 



As one may notice, the projected Hamiltonian contains 4- 
spin and 6- spin interactions in addition to the usual quadratic 
2-spin interactions if it is expressed in terms of the effective 
spin moment ji. One can view these multiple spin terms as 
the interaction between magnetic multipoles (quadrupole and 
octupole) at different sites. Such multipolar Hamiltonians are 
much less familiar than the usual quadratic exchange forms, 
and some caution should be used. In particular, experience 
with similar models shows that such interactions can magnify 
quantum effects, for instance leading to the appearance of a 
quadrupolar phase in the biquadratic case^^. Hence, the naive 
classical approximation - replacing j's by classical vectors - 
is inadvisable, and we will proceed differently below. 



C. Symmetry properties of the Hamiltonian 

Before we move on to discuss the ground state of the 
Hamiltonian H in Eq. (27), we need to have some under- 
standing about its symmetry properties. We start from the NN 
AFM exchange interaction Hex-i- The latter has an apparent 
cubic space group symmetry. The total angular momentum 
J = ^^ji is not conserved, [Hqx-IjJ] 7^ 0. Nevertheless, 
Hex-i surprisingly has a "hidden" SU(2) symmetry. The three 
generators of this global continuous symmetry are defined as 
follows. 



-? 



(28) 



with n = x,y, z. One can readily check that these generators 
commute with Hex-i, 



[G'',Hex-l] =0, 

and satisfy the SU(2) algebra, 



(29) 



(30) 



In addition, the Casimir operator also commutes with 
"Hex-i- The physical meaning of these generators is easy to 
see if one expresses G^'^'^ in matrix form. For a single site. 



Gy = 



Gf = 



1 



1 



(-<t")i4 © (-<t")23 (31) 



= 2('^^)i4®(-^^)23 (32) 



= 2(-'^')i4®(^')23, (33) 
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in which the empty matrix entries are zero and we have ex- 
pressed these generators as the direct sum of two Pauli matri- 
ces, one (0-14) for the subspace of jf = ±3/2 states and the 
other (tT23) for the subspace of jf = ±1/2 states. One in- 
tuitive way to think about these SU(2) generators is that they 
transform the spin components in the jf = ±3/2 subspace to- 
gether with jf = ±1/2 subspace. This is a global symmetry 

ofHex-l. 

Now we consider the other two interactions, Hex-2 and 
^quad- We find that the electric quadrupole-quadrupole inter- 
action Hquad also commutcs with G. On the other hand, the 
ferromagnetic exchange interaction Hex-2 breaks this SU(2) 
symmetry; thus 



[H,G] (X J'. 



(34) 



For J' <C J, y, we have an approximate continuous symme- 
try. 



III. MEAN-FIELD GROUND STATES 



In this section, we study the ground state of the model 



III A 



by get- 



Hamiltonian 1-L in Eq. ( [27] ). We begin in Sec. 
ting some intuition from considering a perturbed model with 
strong easy-axis or easy-plane anisotropy. This starting point 
also has experimental motivation as several ordered double 
perovskites in Tab.|ll]develop such anisotropics that are driven 
by lattice distortions. Armed with the understanding of the 
anisotropic cases, we proceed to analyze the case of cubic 
symmetry using mean-field theory in Sec. IIIB| We consider 



lie 



briefly intermediate strength anisotropy in Sec.[ 

In general. Curie- Weiss mean-field theory consists of de- 
coupling all inter- site interactions to obtain self-consistent 
single-site Hamiltonians. At zero temperature, this is equiva- 
lent to assuming a product form for the wavef unction, i.e. 



(35) 



where the product is over sites, and is an arbitrary j = 
3/2 ket. One calculates the mean-field ground state energy 
as the expectation value of the Hamiltonian in this state, and 
minimizes it. Thus the mean-field approximation can also be 
considered as a simple variational one. 



A. The case with strong anisotropy 



In this subsection, we add to the Hamiltonian 1-L in Eq. ([27 ) 
a strong anisotropic term. 



(36) 



in which, D can be positive or negative, representing easy-axis 
or easy-plane anisotropy, respectively. Although this interac- 
tion is anisotropic in spin space, it still respects the "hidden" 
global SU(2) symmetry. 



1. Easy -axis anisotropy 

Let us start with easy-axis anisotropy D > {). Assuming the 
anisotropy is very strong D ^ J' and D <^ X, which fa- 
vors jf = ±3/2 states, we can safely project the Hamiltonian 
V, into the latter two-dimensional subspace. We then obtain 
the effective Hamiltonian 



E 

(ii> exz.YZ 
J 
4 



J 

4 + 
2 



2 



J'TlTf 



IIV 
~V1 



(37) 



where we have introduced pseudospin-1/2 operators Tj acting 
on the basis j'f = ±3/2 with T? = ±1/2 corresponding to 
j? = ±3/2, respectively. And M is the total number of sites. 
In the reduced space, the original spin vector reduces to the 
pseudospin in the following way: 



(j^j^J")^3(o,o,T^ 



(38) 



Notice that after this projection the interaction on the horizon- 
tal bonds in XY planes disappear in the effective Hamiltonian 
Heff-i. This can be understood in terms of the original orbital 
picture as the easy-axis anisotropy lifts the degeneracy of t2g 
triplets, favoring xz and yz orbitals to be occupied. As a re- 
sult, the above effective Hamiltonian is operating on a bond- 
depleted fee lattice, which is in fact a unfrustrated bipartite 
bcc lattice. 

In Heff-i because of the in-plane anisotropy introduced by 
the FM exchange, the ground state of Heff-i is "antiferro- 
magnetically" ordered in the (T^, T^) plane with an ordering 
wavevector Q = 27r(001). We denote this as the AFM state. 
The corresponding mean-field ground state is just the direct 
product. 



i*w)=ni^'W)' 



where 



V2 



(39) 



(40) 



with an arbitrary phase (j). The arbitrariness of the phase 
comes from the U(l) symmetry of the projected effective 
Hamiltonian Eq. ([37]). However, as discussed in the previ- 
ous section, the continuous symmetry in the original Hamilto- 
nian is broken completely when J' ^ 0. Therefore, the U(l) 
symmetry of Hamiltonian in Eq. ( [37] ) is a by-product of the 
projection. Because we are in the subspace of = ±3/2, the 
orbital occupation is 



(41) 



{{rf^i,yz), {rii^xz), {^i,xy)) = (1/2, 1/2, 0). 



It is also important to note that the ground state in Eq. ( [40] ) 
is not a conventional Neel state as it has a vanishing static 
magnetic dipole moment. 



(^Ij^I^) = 0! 



(42) 
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The 4> dependence only shows up in the spin operators of a 
specific orbital, 

{^\Si,y,\<i>) = -i(-)2-*(cos</.,sin</.,0), (43) 

= ^(-)2-'(cos</.,sin</.,0), (44) 



= (0,0,0). 



2. Easy-plane anisotropy 



(45) 



with 



V2 



■\)] 



(49) 



This is once again an antiferromagnetic state, and to distin- 
guish it from the one which obtains for Ising anisotropy, we 
denote it AFM'. The defining difference of the AFM' and the 
AFM state discussed previously is that, the former has a non- 
zero dipole moment, while, at least within mean field theory, 
the latter does not. 



Now we consider easy-plane anisotropy D < {). We also 
assume the anisotropy is very strong \D\ ^ J^J'^V and 
\D\ <C A, which favors jf = ±1/2 states and obtain the ef- 
fective Hamiltonian after projection into the jf = ±1/2 sub- 
space. 



^eff-2= E 9('^^' 



{ij)exY 
1 



E 

{ij)exz 



J_ 
36" 



T, ■ T, 



_rpXrpX 

6^ ^ 



rpyrpy 
18 ' ^ 



_rTiZrTiZ \ 

2 i .7 / 



E 

(ij)GYZ 
1 

6 



36' 



rpZrpZ \ 

i 3 ) 



J 

4 



3J^ 
2 



nv 



6 ' ^ 



). 



(46) 



Here the pseudospin-1/2 operator T is acting on the subspace 
of f = ±1/2 with = ±1/2 representing = ±1/2, re- 
spectively. In the reduced spin space, the original spin vector 
is reduced to the pseudospin in the following way: 



if 



(2T^,2T^T^ 



(47) 



B. The cubic case 

Having understood the cases with strong easy-axis and 
easy-plane anisotropics, let us now turn to the Hamiltonian 
H with cubic lattice symmetry in Eq. ( [27] ). Both ground 
states of the Hamiltonian with strong easy-axis or easy-plane 
anisotropy comprise two-sublattice structure with an ordering 
wavevector equivalent to Q = 27r(001). It is therefore natu- 
ral to guess that the same two sublattice structure is also ob- 
tained in the cubic case. While we have not proven this, we 
have investigated more general mean-field ground states, al- 
lowing for much larger unit cells, but found in every case that 
the minimum energy is found for the two-sublattice configura- 
tion. Therefore, in what follows, we assume the two sublattice 
structure with ordering Q = 27r(001) (which is equivalent to 
27r(100) and 27r(010) in the cubic case). We make no fur- 
ther assumptions, and minimize the energy with respect to an 
arbitrary wavef unction on each of the two sublattices. The 
resulting variational phase diagram is depicted in Fig.[T] 



1. Antiferromagnetic (AFM) state 



We can now find the mean-field ground state of this Hamil- 
tonian. For an effective S = 1/2 model of this type, this 
is equivalent to the classical approximation. Classically, we 
can find the minimum energy states by the Luttinger-Tisza 
method. This amounts to looking for the eigenvalues of the 
spin Hamiltonian regarded as a quadratic form, and seeking a 
classical spin solution which is built of a superposition only 
of those eigenvectors which have minimum energy eigenval- 
ues. The result in this case is that, for < J' < J, there are 
two classes of solution, all collinear spin states. These are: 
(i) states with Q = 27r(100) and the pseudo-spin axis in the 
yz plane, and (ii) states with Q = 27r(010) and the pseudo- 
spin axis in the xz plane. As for the easy-axis case, there is 
an accidental degeneracy of spin orientations within the plane 
normal to Q. Note that while the pseudospin orients freely 
along a circle in this plane, the magnetization orients along 
an ellipse due to the factor of 2 in Eq. ( [47] ). One readily ex- 
presses the ground state in the unprojected Hilbert space. For 
example, taking Q = 27r(010) and pseudospin pointing along 
X direction, then 



i*>=ni^*' 



(48) 



In Fig.[T] for small J' /J and V/J, we find a phase, denoted 
AFM, which is the natural continuation of the AFM phases 
encountered in the anisotropic limits. Here, as in those cases, 
the states on the two sublattices are simply related by a time 
reversal transformation, and indeed the ground state has the 
same form as that found in the easy-axis case, Eqs. ([39|),(|4Q|). 
The appearance of time-reversed pairs of sites is natural, since 
the largest interaction, Hex-i , is dominated by the spin-spin ex- 
change term. Interestingly, one finds that the ground state has 
a continuous degeneracy: the phase (j) in Eq. ( |4Q| can be ar- 
bitrary. Since the Hamiltonian with non- vanishing J' has no 
continuous symmetry, this degeneracy appears to be acciden- 
tal. Since it has the same form as we found in Sec. |III A 1 1 
we continue to use the label AFM here for this state (and in 
Fig.[l). 



2. Ferromagnetic 110 ( FMl 10) state 

With large J^J and V/J, the orbital-orbital interaction 
has more weight in the Hamiltonian and the nature of the 
ground state changes. One should note that even the pure 
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orbital-orbital interaction is not trivial and classical, since the 
orbital occupation numbers no longer commute after projec- 
tion down to the j = 3/2 quadruplets. However, in mean- 
field theory one may still treat the expectation values classi- 
cally. Note that the largest terms in the orbital-orbital interac- 
tion are those which are diagonal in the orbital basis, namely, 
the second term in Eq. ([19]) and the second term in Eq. pT] ). 
To minimize the diagonal orbital interaction like fii^xz^j.xz, a 
schematic recipe is to maximize fii^xz while minimizing fij^xz- 
This is necessary because one cannot minimize both fii^xz and 
fij^xz simultaneously, since, due to the single-occupancy con- 
straint, the other diagonal terms such as fii^yzfij^yz would then 



be increased. Since the occupation numbers of the same or- 
bital must be taken different on different sites, and these oc- 
cupation numbers are time-reversal invariant, the states on the 
two sublattices cannot be time-reversed counterparts. Con- 
sequently, there is a competition between the orbital-orbital 
interactions (Hex-2 and Hquad) and the nearest neighbor anti- 
ferromagnetic exchange interaction (Hex-i). In the large J' / J 
and V/ J region, when the orbital interactions dominate, we 
find however that time-reversal symmetry is still broken, and 
since these states are not composed of time-reversed pairs, the 
result is an uncompensated net ferromagnetic moment. 

In the majority of phase space, we find the ground state is 
characterized by three parameters, r, and 



|^^>FMUO = ^(e^^Mj^ = i> + e^(^-^^)b- = -\)) + \^{e^'^^\f = §> + \f = -§» (50) 



|^b)fmuo = -^(-e-^M/ = i) +ie*(^-^^)|r = -i)) + ^ ^{-ie-^f^lf = |) + \f = -§)) , (51) 

in which, "A" and "B" represent the two sublattices, and r, (/)i and (/)2 are determined by minimizing the mean field energy. Note 
that in Eq. ( [5T] ) the three parameters r, (/)i and (/)2 are uniquely determined by J' / J and V/ J. So the orbital occupations can be 
readily generated, 

1 Ty/ 1 — 1 vy/ 1 — 

{{riA.yz), {riA.xz), {nA,xy))vMnQ ^ ~ y 71 — cos^i, 2 ~ y ^ — COS 01, -^), (52) 

1 T\/ 1 — 1 1 — 2r^ 

{{riB.yz), {riB.xz), {riB,xy))vMnQ = (2 ~ y ^ 71 — cos 01, 2 ~ y — cos 01, — ). (53) 

It is interesting to see the spin vectors of two sublattices are symmetric about [110] direction, 

0'a)fmiio = r (^\/3 - 3r2 cos (0i - ^2) +rcos(20i - 202), \/3 - Sr^ sin (0i - 02) - rsin(20i - 202), 0^ (54) 

0'b)fmiio = r (^-\/3 - 3r2 sin (0i - 02) + rsin(20i - 202), -\/3 - 3r2 cos (0i - 02) - rcos(20i - 202), 0^ (55) 

so the system has a non- vanishing net spin polarization, that is 

T r / 1 

\{3a^ jEhmno ^ 2 i 3r2(cos(0i - 02) - sin (0i - 02)) +r(cos(20i - 02) + sin (20i - 02))J (1,-1,0) . (56) 

This direction of polarization is equivalent to [110] by a 90 degree rotation, so we denote this a FMllO state. It occupies the 
corresponding region in Fig.[T] 

3. Ferromagnetic 100 (FMIOO) state 

Between the AFM and FMl 10 states, a narrow region of intermediate phase intervenes (see Fig.[T]). Numerically we find this 
phase is characterized by two parameters ri and r2 , 



|^^>FMioo = ^(1/ = \) + \f = -i» + {\f = 1) + \f = -§)) (57) 



|^b)fmioo = ^i-\f = h) + \f = -!)) + \J^(-\f = !) + \f = -§))• (58) 

The parameters ri and r2 are determined by J^J and V/J and in this intermediate phase ri ^ r2. A second ground state is 
obtained then by interchanging ri and r2 . The orbital occupation numbers and spin vectors are give by 



V^V i'2 3 ' 1' 3 

{{riB.yz)^ {nB,xz): {nB,xy))FMlOO = ( 2 ~ ^ ^ V ^ ~ ^2 ^ 2 ~ ^ ~ V ^ ~ ^2 ^ "3^) ^^^^ 
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and 



FMIOO 



FMIOO 



^ +ri./3-3rf,0,0 



r2 



3-3r2 0,0 



(61) 



(62) 



We see that the net spin polarization is along the [100] direction. Due to cubic symmetry, all possible [100] directions are 
possible. By analogy with the previous phase, we denote this phase FMIOO. It occupies the narrow region shown in Fig.[T] 



4. Transitions 

The intermediate FMIOO state in Eq. ( [58] ) is smoothly con- 
nected to the AFM state, which is obtained by setting ri = 
r2 = 0. By contrast, it cannot be connected to the FMllO 
state. This indicates that the transition between FMIOO 
to AFM is continuous while the transition from FMIOO to 
FMllO is first-order. Indeed, this can also be clearly seen 
from the behavior of the spin and orbital order parameters 
across these transitions (see Fig.[3|. Both the spin and orbital 
order parameters jump when the system goes from FMIOO 
phase to FMl 10 phase. 

The approach to the line J' = 0, while not a transition per 
se, does represent a change of behavior. In particular, on this 
line, the hidden SU(2) symmetry is restored, and new ground 
states may be obtained from the above three phases by SU(2) 
rotations. This allows, for instance, for AFM states to develop 
with non- vanishing magnetic dipole order in the ground state. 
However, for arbitrarily small J' the SU(2)-induced degener- 
acy is broken and the results quoted above hold. 




0.23 0.24 



0.25 0.26 




0.25 0.26 



C. Intermediate anisotropy 

We briefly address here the evolution of the ground 
states with increasing \D\ between the cubic and strongly 
anisotropic limits. For easy-axis anisotropy, D > 0, this 
process is relatively simple. The AFM phase (with Q = 
27r(001)) is favored by this sign of anisotropy, and therefore, 
with increasing D, it expands at the expense of the FMllO 
and FMIOO states. Indeed, for very large D, ferromagnetic 
states occur only for unphysically large J' and V. 

In the case of easy-plane anisotropy, D < 0, the phase dia- 
gram is more subtle. For weak \D\, the main effect is to break 
the symmetry between the formerly equivalent [100] wave vec- 
tors. In this case, states with minimal (j^) are favored, which 
prefers Q = 27r(100), 27r(010) rather than Q = 27r(001). 

If we begin in the AFM state for D = 0, this aligns 
the pseudospin in the plane normal to this wavevector. The 
phase degeneracy which obtains for the cubic case is bro- 
ken by the anisotropy, and a definite alignment is obtained. 
Moreover, as states with = ±1/2 are increasingly mixed 
into the ground state, a non-vanishing dipole moment, pro- 
portional to the pseudospin, is induced. The magnitude of 
this staggered magnetization grows continuously with \D\, 
eventually sls D — oo, approaching the value obtained in 
Sec. [Ill A 2[ This local moment is oriented in the plane nor- 







as***" 


V 


' 0A+iB)/2 




V 


' Oa+;b)/2 







0.21 0.22 0.23 0.24 0.25 0.26 



V 

FIG. 3. (Color online) Upper graph: fixy versus V; middle graph: 
fixz and fiyz versus V; lower graph: the net spin polarization per site 
versus V. In all three graphs, J = 1 and J' = 0.2, and T = 0. 
The assignment of A and B sublattices is arbitrary. Here we take the 
choice given in the text. 



mal to Q, and can take values distributed over an ellipse in 
this plane. In the large \D\ limit, this ratio of the major (per- 
pendicular to Q and to z) and minor (z) axes of the ellipse 
approaches 2, corresponding to the accidental degeneracy dis- 
cussed in Sec. IIII A 21 . Because the state for non-zero D 
evolves smoothly into this limit, and has a non-zero local mo- 
ment, we denote it an AFM' state, following the earlier nota- 
tion. 
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Beginning in the FMllO state at = 0, one observes two 
subsequent transitions. First, small \D\ orients the magnetiza- 
tion normal to Q. For concreteness consider Q = 27r(100), in 
which case one obtains a ferromagnetic magnetization of the 
form m = (0, mi, 7712). For D = 0~, mi = m2, but sub- 
sequently m2 decreases such that m2 < mi. We denote this 
state FMllO*. An example for the orientation of net polar- 
ization in the FMl 10* by varying V is given in Fig.|4] Even- 
tually once some critical anisotropy is reached, m2 vanishes 
continuously. At this point the magnetization is aligned along 
the (010) axis. For yet larger anisotropy, eventually the fer- 
romagnetic magnetization vanishes entirely, and the ground 
state switches to the AFM' state. An example of the transi- 
tions from FMl 10* to FMIOO then to AFM' by varying the 
easy-plane anisotropy D is given in Fig.|5] 

Finally, starting in the FMIOO state, the magnetization im- 
mediately switches to the (010) direction. This is the same 
phase as the intermediate phase observed starting from the 
FMl 10 phase. Thus with further increase in anisotropy, the 
ground state switches to the AFM' state. 

One may also visualize the evolution of the ground states 
with anisotropy by considering planar phase diagrams at fixed 
D. With increasing positive D (Ising anisotropy), the AFM 
state is stabilized, and simply expands in the — V plane, 
pushing the FMIOO and FMl 10 states outward. For increas- 
ing negative D, apart from the fact that the AFM and FMl 10 
states evolve into the AFM' and FMl 10* states, the behav- 
ior is similar: the AFM' state expands at the expense of the 
ferromagnetic states. 




V 



FIG. 4. (Color online) The orientation of the net spin polarization 
for the FMl 10* phase, at T = 0. ^ is angle between the net spin 
polarization and nearest [110] direction. In the figure, J' = 0.2, D = 
—0.05, J = 1. cos increases from 1 / >/2 for the FMIOO phase to 1 
for the FMl 10* phase as V goes through the phase transition point. 



IV. MULTIPOLAR ORDERS AND T > BEHAVIOR 

A. Order parameters 

In this section, we extend the analysis of the previous sec- 
tion to non-zero temperature. To do so, we employ the usual 
extension of mean field theory to include thermal fluctuations. 
To characterize the phases encountered in this treatment, it is 




D 



FIG. 5. (Color online) The net polarization versus the easy plane 
anisotropy, at T = 0. The upper curve (in red) is for the y component 
of the net polarization. And the lower curve (in blue) is for the z 
component of the net polarization. The ordering wavevector is Q = 
27r(100). In the graph, f = V = 0.4 and J = 1. When < |D| < 
0.36, the system is in FMl 10* phase; when 0.36 < \D\ < 0.41, the 
system is in the FMIOO phase; when \D\ > 0.41, the system is in 
AFM' phase. 



natural to introduce several types of order parameter. First, on 
a single site i, we may measure the dipole moment, which is 
proportional to ji. However, we may also measure the next 
two multipoles: the quadrupole moment, proportional to 



QT = OTjT) 

and the octupole moment 



(63) 



(64) 



A typical magnetic state has a non-vanishing local dipole mo- 
ment, which inevitably induces some higher multipole order 
parameters (see below). However, one sometimes encounters 
purely multipole states, in which (jf ) = but Qf^ and/or 
Of ^ are/is non- vanishing. 

The components of these tensor can be decomposed into 
irreducible representations of the cubic group (characterizing 
the symmetry of the ideal double perovskite structure). This 
decomposition is described fully in Table |T| Here we note in 
particular the two-dimensional Fs representation 



= m' 



1)) 



(65) 



which are analogous to the Cg orbitals in atomic physics. The 
remaining three independent components of Qf ^ (jfji-^jiJi 
etc.) form a three-dimensional representation analogous to the 
t2g orbitals, but do not appear in our analysis. 

Another important way to break up the tensor order param- 
eters is into combinations which appear in the spin Hamilto- 
nian. Specifically, these are the orbital occupation operators, 
^i,yzj^i,xz^^i,xy, and the orbitally-resolved spin operators, 
^i,yz^ ^i,xz^ ^i,xy Thcsc Can be expressed in terms of the 



11 



Moment Symmetry 



Dipole 



Quadrupole 



Octupole 



r2 
r4 



Operator 



My = f 



= [3(r)'-j']/V3 

Qyz ^ 



tI = iff - [jyifr + myy]/2 

z\3 



Ta = U 



\j-{j-y + {jy)^j ^]/2 



TABLE L Multipole moments within a cubic Fs quartet. Bars over 
symbols indicate the sum with respect to all the possible permuta- 
tions of the indices, e.g. 'j^W? = f{fY + fff + {fff- 
Adapted from Ref.Hand Ref.31. 



1 2 2 

-y 



multipoles describe above. For the occupation numbers, 
1 1 

""^'"^ " 3 ^ 

_ 1 1 
""^'^^ " 3 + 

_ 1 1 
""^'^^ " 3 " 



1 2 2 

J- -y 



(66) 



The orbitally -resolved spins decompose as 



i,yz 



SI 

QX 

qy 



qX 



dy 

i,xy 



St 



tjXy 



1 

15^ 
2 

15 
2 

15 
2 

15 
1 

15 
2 

15 
2 

15 
2 

15 
1 

15-^ 



15 ''^ 
{y _i_ —T^ -I - 



j7 + 



j7 + 



15 
15 
15 
15 
15 
15 
15 



3V15 
1 



3V15 
1 



3V15 



1 



3V15 
1 



3V15 
1 



3V15 



rpy 



(67) 



B. Cubic system: phases 

We first discuss the phases occurring in the cubic system 
at T > 0. The ground states discussed earlier are all stable 



to small thermal fluctuations, and hence persist at low tem- 
perature. Thus we expect, broadly speaking, an antiferro- 
magnetic (AFM) and ferromagnetic (FMllO/FMlOO) region 
at low temperature. Of course, at temperatures much larger 
than J, one has a disordered paramagnetic phase. Interest- 
ingly, an additional phase appears at intermediate tempera- 
ture. This is a non-magnetic quadrupolar ordered phase. 

To see how this arises, we describe the mean-field pro- 
cedure and its results. Mean field theory is formulated in 
the usual way. We self-consistently decouple interactions be- 
tween different sites i and j as follows: 



(68) 



where Oi and Oj are two operators at site i and j, respectively. 
These operators are nothing but the orbital occupation num- 
bers and orbitally resolved spins, which are related to the mul- 
tipolar operators by Eqs. ([66])- ( [67] ). Decoupling all pairwise 
interactions between sites in this way, we then obtain a set of 
single-site problems for each j = 3/2. Note that these single- 
site problems involve not just the usual Weiss exchange field, 
but also "multipolar fields", which act as effective second and 
third order spin anisotropics. The mean-field equations deter- 
mine self-consistent values of the orbital occupation numbers 
and orbitally resolved fields. As it is straightforward to for- 
mulate the mean-field equations, and solve them numerically, 
we do not give the details of these calculations here. 

A distinct class of solutions describes each phase. For the 
antiferromagnetic phase, we find the following operators are 
non-zero: 



{ji) = ±n{ui,U2,0), 



(69) 
(70) 
(71) 
(72) 



where we have taken Q = 27r(0, 0, 1) for concreteness, and 
the upper and lower signs refer to the A and B sublattices, 
respectively. The parameters n, g',^^, and t/3 are positive at 
all T > in the AFM phase. However, note that n vanishes 
in the limit T 0, in agreement with the vanishing dipole 
moment discussed earlier for the AFM ground state. 
In the FMl 10 state, the non-zero expectation values are: 

(j,) =m(l,l,0)±n(l,-l,0), (73) 

{Qf)=q. (74) 

{Qf-'')=Tq\ (75) 
(T,,,) =t«(l,l,0)±t«(l,-l,0), (76) 
(T,,^) =t^(l,-l,0)±t^(l,l,0), (77) 

where again we took Q = 27r(0,0, 1) and the upper/lower 
signs refer to the A/B sublattices. In this case the parame- 
ters m^n^q^q' ^ta^ia^tf3^if3 are all non-zero at temperatures 
within the FMl 10 phase including T = 0. 

The third ordered phase dominating the phase diagram is 
the quadrupolar one (For the purposes of this section, we ig- 
nore the FMl 00 phase, which extends into a narrow region 
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of ferromagnetic state with variable polarization direction at 
T > 0, as it occupies a very small volume of the phase di- 
agram). In the quadrupolar state, there is only a single non- 
vanishing order parameter: 



{Q- ) 



(78) 



Let us discuss the symmetries of these three states. In the 
AFM and FMllO phases, time reversal symmetry is broken. 
However, the net magnetization vanishes in the AFM state. 
In the AFM state, this is guaranteed by invariance under the 
combined operations of translation (such as by (0, 1/2, 1/2), 
which interchanges the A and B sublattices) and time-reversal. 
No such symmetry can be combined with time-reversal in the 
FMl 10 case. Various point group symmetries are also present 
in the AFM and FMl 10 phases, but we do not describe this in 
detail. 

In the quadrupolar case, time-reversal symmetry is unbro- 
ken, which is sufficient to require the dipolar and octupolar 
order parameters to vanish. Only point group symmetries are 
broken by the quadrupolar order. Four- fold (C4) rotations 
about the xory axes, and three-fold (C3) rotations about [111] 
axes are broken in this state. While the C4 rotation about the 
z axis is also broken, the combination of this C4 rotation and 
a translation exchanging the A and B sublattices remains a 
symmetry of the quadrupolar state. 

A standard classification scheme for quadrupolar states is 
to examine the examine the eigenvalues of the Qf^ ma- 
trix. These must sum to zero because the matrix is trace- 
less. States in which there are only two distinct eigenvalues, 
i.e. eigs((5) = {q^Qj—'^q} are called nematics, and corre- 
spond to the situation in which one principal axis is distin- 
guished from the other two, which remain identical. In the 
most general case, there are three distinct eigenvalues, i.e. 
eigs{Q) = {gi, ^2, -Qi - ^2}, with qi q^. This is called a 
biaxial nematic, and is a state in which all three principal axes 
are distinct. The quadrupolar state obtained here is such a bi- 
axial nematic. Physically, the local susceptibility in this state 
takes distinct values Xkfcai^ Xkfcai' Xkfcai ^^1^^ along each 
of the axes. However, note from Eq. ( |78| ) that the quadrupo- 
lar order parameter changes sign between the two sublattices. 
Thus we should properly call this state an antiferro-biaxial ne- 
matic. Due to the staggered ordering, the bulk susceptibility 
does not distinguish all three axes. Instead, there are only 



two distinct components, 



X^^ 7^ X^^- The difference 



between the two components of the susceptibility serves as a 
simple macroscopic means to observe quadrupolar ordering. 



C. Cubic system: phase diagram and transitions 

By solving the mean-field equations numerically, we have 
determined the phase diagram for the cubic case. Parts of it 
can be understood analytically. Suppose that the transitions 
from the high temperature normal phase to the quadrupolar 
and AFM phases are second order. This appears to be always 
true for the quadrupolar phase, while it true for the AFM for 
most parameters, but weakly violated in some regions. With 



this assumption, we can determine the critical temperatures 
for these transitions by the usual condition of marginal stabil- 
ity (vanishing of the quadratic term in the Landau theory) of 
the free energy. We find that the critical temperature for the 
quadrupolar state is 



43V + 18 J' - 3 J 



18 



(79) 



and that for the AFM state is 



^AFM J + 10 J' + x/73 J2 + 164 J J' + lOOjjr 

= . (80) 

Without the assumption that the transitions are continuous, the 
critical temperature could be higher. Thus Eqs. ( 79|80| ) give 
lower bounds for the transition temperatures, strictly speak- 
ing. Extending the two-dimensional T = phase diagram in 
Fig. [T] into a third dimension of temperature, the quadrupolar 
phase occurs "above" the portion for which T^^^^ —^^^^ 
which implies V > with 



Vr. 



7 J - 26 J' + V'73J2 + 164JJ' + 100( J' 
86 



f\2 



(81) 



The curve Vc{J') defines an almost straight line in the 2d 
phase diagram, as shown in Fig. [6] In fact, Eq. ( [ST] ) slightly 
underestimates as it assumes the normal to AFM transi- 
tion is continuous, when it is in fact weakly first order in this 
vicinity. However, the true Vc found numerically is only a few 
percent larger. For V < no quadrupolar phase occurs. 
Instead, the first ordering transition from high temperature is 
into the AFM state. This is true even when the ground state is 
ferromagnetic, so that in this case (when V < Vc) the system 
first orders into the AFM state, and then at lower temperature 
switches to the FMl 10 phase. 

We now discuss the nature of the transitions. The quadrupo- 
lar ordering transition is, as already mentioned, continuous 
(see Figs. [8] |9]). It is described by a single scalar order pa- 
rameter, (equal to q' in Eq. (|78|)), for each of the three [100] 
wavevectors, describing the associated staggered quadrupole 
moment. Formally, 



(82) 



where x\ = = = Zi. According to symmetry, 

the Landau free energy for (j)a has the same form as that for an 
0(3) magnetic transition with cubic anisotropy. Beyond mean 
field theory, this is believed to support a three dimensional 
transition in the 0(3) universality class. 

The transition from the normal to the AFM state is contin- 
uous in mean field theory for small V (see Fig. [7]), becoming 
weakly first order for larger V, close to Vc where the interme- 
diate quadrupolar phase emerges. The normal- AFM transition 
is characterized, for each of the three wavevectors, by a two- 
component primary order parameter, which could be taken to 
be the two components of T^^q, normal to Q. In principle, 
the degeneracy of the ordering pattern within this "XY" plane 
normal to Q is, as we have remarked, accidental, and should 
be removed by additional effects. We do not, however, ob- 
serve this degeneracy lifting within mean field theory for the 
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FIG. 6. (Color online) Zero temperature two-dimensional phase di- 
agram in the cubic case (same as Fig. [TJ, overlaid with the regions 
of different T > behavior. Dashed lines are obtained from mean 
field numerics. The solid line (in orange) is defined by Eq. jsT] ). 
In region I, there is a single transition to the AFM state. In region 

II, the system supports an intermediate temperature quadrupolar or- 
dered phase. The transition from the normal state to the quadrupolar 
phase is second order. In this region, there is a first order transition 
to the AFM phase on further cooling, for parameters such that the 
latter is the ground state. Otherwise, the lower temperature transi- 
tion is to a ferromagnetic (predominantly FMllO) state. In region 

III, the system first turns from normal phase to the AFM state then to 
a ferromagnetic (predominantly FMllO) state on further cooling. 




FIG. 7. (Color online) Temperature dependence of order parame- 
ters fox J' / J = 0.2, V/ J = 0.1. For these parameters, there is 
a direct, continuous, normal to AFM transition, dX T / J ^ 0.38. 
The three curves show: squares (red online) | (Ta,^ — Tb,^) |/2, cir- 
cles (blue online) (Qjf + Q^b )/2, and diamonds (yellow online) 
IOa — 3b)\I'^- Note: in Figs. 7-10, the symbols are not the data 
points (which are much more dense) - they simply label the different 
curves. 



present model. With the degeneracy, the transition should be 
therefore described by the free energy for some six component 
order parameter. As we do not understand the degeneracy lift- 
ing mechanism at present, we do not attempt here to establish 
the true critical properties for this transition (when it is con- 
tinuous) with fluctuations taken into account. 

The quadrupolar to FMl 10 transition is continuous in mean 
field theory (see Fig.|9]). This could be anticipated by exam- 
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FIG. 8. (Color online) Temperature dependence of order parameters 
for J' I J = 0.1, y/ J = 0.3. For these parameters, there is a contin- 
uous normal to quadrupolar transition, dXT/ J ^ 0.65, followed by 
a first order transition to the AFM state at T/J ^ 0.37. The four or- 
der parameters plotted are: squares (red online) | (Ta,^ — TB,ct) |/2, 
circles (blue online) (Q^f + Q^b )/2, diamonds (yellow online) 
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FIG. 9. (Color online) Temperature dependence of order parame- 
ters for J' I J — 0.3, V I J — 0.3. For these parameters, there is a 
continuous normal to quadrupolar transition, diiT/ J ^ 0.85, fol- 
lowed by a continuous quadrupolar to FMllO transition dXT/ J ^ 
0.33. The four order parameters plotted are: squares (red online) 
|(TA,a - Ts,a) 1/2, circles (blue online) {Q^/ + Qs')/2, dia- 
monds (yellow online) |(ja + jfs)|/2, and triangles (green online) 



ining the form of the order parameters in the FMllO phase. 
We note that the antiferro-biaxial nematic order parameter of 
the quadrupolar state is already non- vanishing in the FMllO 
phase. Hence, we might naturally expect, upon heating, 
that thermal fluctuations first restore time-reversal symmetry, 
yielding the quadrupolar phase, before fully restoring all sym- 
metry in the normal state. To determine the nature of the or- 
der parameter for this transition, note that the wavevector and 
local anisotropy axes are already established in the quadrupo- 
lar state. Hence the direction of the uniform and staggered 
magnetizations are already determined, up to a sign and inter- 
change, above the transition. For instance, for the quadrupolar 
state in Eq. ([78]), with Q = 27r(001), the uniform magneti- 
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Qf field: 
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T/J 

FIG. 10. (Color online) Temperature dependence of order parame- 
ters for J' /J = 0.40, V I J — 0.05. For these parameters, there is a 
continuous normal to AFM transition, dXT/ J 0.49, followed by a 
continuous AFM to FMllO transition at T/J ^ 0.34. The four or- 
der parameters plotted are: squares (red online) | (Ta,^ — Tb,ck) |/2, 
circles (blue online) (Qjf + Q^^ )/2, diamonds (yellow online) 
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zation can lie along ±(110) and the staggered magnetization 
along ±(110), or vice- versa. Thus the symmetry breaking 
from the quadrupolar to the FMl 10 state is described by two 
Ising order parameters. We therefore expect this transition, 
beyond mean field theory, to be similar to that of an Ashkin- 
Teller or similar models. 

The quadrupolar to AFM transition appears strongly first 
order (see Fig. [8]). This is in agreement with the expectations 
of Landau theory, as the symmetry of the AFM phase is not a 
subgroup of the symmetry of the quadrupolar one. In terms of 

2 _ 2 

order parameters, this is evident since {Q^ ^ ) is non-zero 
in the quadrupolar phase but zero in the AFM one, while the 
magnetic order parameters are zero in the quadrupolar phase 
but non-zero in the AFM. Fine tuning of the free energy would 
be required to arrange both these types of order to change at 
the same temperature in a continuous fashion. 

In region III one encounters a transition from the AFM to 
FMl 10 state. This appears to be continuous in mean field the- 
ory (see Fig. [TO]). One can understand this by noting that the 
AFM solution can be regarded as a subset of the FMl 10 one, if 
the unit vector (iii, ii2, 0) is taken to be along (1,-1,0). Then 
the transition to the FMl 10 is described by the emergence of 
a non-zero m. Like the normal to AFM transition, because 
we have not understood the degeneracy -breaking mechanism 
in the AFM state, we do not attempt to analyze this transition 
beyond MFT. 



D. Effects of anisotropy 

We now consider the effects of anisotropy on the T > 
phase diagram, focusing on the case of weak \D\ <C J^J'^V. 
We have already considered the effects of D on the AFM and 



FMl 10 states in Sec. |IIIB 4| We saw that easy-axis anisotropy 
favors states with the wavevector Q parallel to the z (easy) 
axis. This is because the anisotropy couples directly to the 
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const. 
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(83) 



Both the AFM and FMl 10 states have a non-zero and constant 
expectation value of Qf ^ , which is maximized in this orien- 
tation. Conversely, easy-plane anisotropy favors states with 
the wavevector Q perpendicular to the z (hard) axis, for the 
same reason. 
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FIG. 1 1 . Temperature dependence of the free energy difference be- 
tween states with wavevector parallel and perpendicular to the Ising 
axis, in the presence of a weak anisotropy = 0.05. Here 
f/J = V/J = 0.3. Solid line: AF = F{Q \\ z;D = 
.05) - F{Q ± z]D = .05). Dotted line: AF = F{Q \\ z;D = 
— .05) — F{Q ± z;D = —.05). One sees that in the quadrupolar 
phase, both signs of anisotropy favor the wavevector aligned with the 
z axis. In the FMl 10 phase, however, this is favored only for D > 0. 
For D < (easy plane anisotropy), the state with wavevector normal 
to z is preferred. Note also that the energy difference is much larger 
in the FMl 10, consistent with the expected linear and quadratic de- 
pendence on D in FMl 10 and quadrupolar phases, respectively. 

We now repeat this analysis for the quadrupolar state. 
Here the situation is more subtle because (Qf^ ) vanishes 
in the quadrupolar state. Moreover, the cubic rotations (e.g. 
(Qf^ )), while not vanishing, give zero net contribution due 
to the opposite signs on the A and B sublattices. This means 
that the splitting of the different wavevector states vanishes 
at linear order in D. There is instead a quadratic contribu- 
tion, which, numerically, we find favors the states with Q 
parallel to z (see Fig. [TT]). Being quadratic in D, this same 
configuration is favored /(9r both the easy-axis and easy-plane 
case. Thus we have the interesting situation that for easy- 
plane anisotropy, the wavevector orients parallel to z in the 
quadrupolar phase, but perpendicular to z in the low temper- 
ature phase. Note that the quadrupolar phase remains distinct 
from the normal phase even with non-zero I), as it continues 
to break symmetries, notably translational invariance. 



E. Magnetic susceptibility 

In this subsection, we discuss the magnetic response at T > 
0, which is an important indicator, especially of the quadrupo- 
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FIG. 12. (Color online) Inverse susceptibility at the normal to AFM 
transition for = V — O.IJ. Blue (upper) curve: 1/xxx, red 
(lower) curve: 1/xzz- 



lar ordering transition. At high temperature, of course, one 
observes Curie- Weiss behavior. For the general Hamiltonian 
with anisotropy D, there are two different Curie- Weiss tem- 
peratures, for fields parallel and perpendicular to z: 
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(84) 



These are obtained from the high temperature expansion of 
the susceptibility up to 0(1/T^). These expressions may be 
useful in extracting exchange constants from experiment. In- 
terestingly, if one calculates the powder average average sus- 
ceptibility, the contributions of the anisotropy cancel at this 
order and the Curie-Weiss temperature measured in this way 
is independent of D. It is also interesting to note that, in the 
region of larger V/ J and small J' / J, one obtains a ferromag- 
netic ground state with an antiferromagnetic (negative) Curie- 
Weiss temperature. 

On lowering temperature, the susceptibility shows distinct 
behaviors in the different parts of the phase diagram. We focus 
here for simplicity on the cubic system, starting with region I. 
Here the susceptibility displays the usual cusp associated with 
antiferromagnetic order, at the normal to AFM transition. The 
inverse susceptibility is plotted in Fig. 12 for J' = F = O.IJ, 
in the midst of region I. It shows a minimum at the transition, 
and pronounced curvature below the transition temperature, 
saturating to a large constant value in the T ^ limit. We 
note that the large zero temperature susceptibility is not re- 
lated to gapless excitations, but is a general consequence of 
strong SOC, and should be expected in all parameter regimes 
of this model. 

Next consider region II. Here, one observes a cusp at the 
normal to quadrupolar transition. This cusp is, however, rather 
different from the one just mentioned. Specifically, it is not a 
minimum of 1/x, and instead separates two distinct "Curie- 
Weiss" regimes in which 1/x is linear but with different, pos- 
itive, slopes (i.e. different effective magnetic moments). The 
presence of a lower temperature Curie- Weiss regime is a sig- 
nature of quadrupolar order. This is because the quadrupo- 
lar mean field splits only the point group degeneracy of the 
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FIG. 13. (Color online) Inverse susceptibility for for J' — 0.2 J, 
V — 0.3 J. Blue (lower) curve: 1/xxx, red (upper) curve: 
yellow (middle) curve: l/xpowder- For these parameters the 
quadrupolar transition is at T/ J 0.75, and the ferromagnetic tran- 
sition is at T/J 0.23. 



spins, but preserves a local Kramer's doublet. This doublet 
gives rise to a Curie law. An example is plotted in Fig. [13 
As the quadrupolar order lowers the symmetry of the system 
to tetragonal, we see actually two different effective moments 
in susceptibility parallel to the wavevector Q (Xzz) and per- 
pendicular to it (Xxx = Xyy)- We observe that the effective 
magnetic moment seen in Xxx is typically enhanced in the 
quadrupolar phase, while it is suppressed in Xzz^ both rela- 
tive to the isotropic effective magnetic moment in the normal 
phase. 

At still lower temperature, one encounters the ferromag- 
netic phases. Here of course the susceptibility for the easy 
directions diverges. Focusing on the dominant FMllO phase, 
one sees that since the easy direction is in the (001) plane se- 
lected by the quadrupolar order, Xzz does not diverge, but Xxx 
and Xyy do. 



V. 



BEYOND MEAN-FIELD: SPIN WAVES AND 
NON-MAGNETIC GROUND STATES 



In Sec. VA we obtained the mean field phase diagram. 
Here we consider quantum effects beyond mean field. We 
first consider spin wave fluctuations, and obtain the collec- 
tive mode spectrum in linear spin wave theory. From this, 
we obtain the quantum correction to the order parameter, and, 
in the ideal case of J' = V = 0, we will see that this is 
very large and invalidates the mean field theory in the vicin- 
ity of this parameter regime. This suggests the possibility of 
very different states dominated by quantum fluctuations. We 
then explore this possibility, considering some candidate non- 
magnetic ground states of our model. 

First, we consider the quantum ground states of pairs of 
sites, unveiling a pseudo-singlet structure, analogous to the 
5* = singlet ground states for pairs of antiferromagnetically 
interacting spins with SU(2) symmetry. This leads naturally 
to the possibility of "valence bond" states built from these 
pseudo-singlets. We consider both a static. Valence Bond 
Solid (VBS) state, and states in which the valence bonds are 
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fluctuating, in which case we obtain a Quantum Spin Liquid 
(QSL) state. 

Very Httle theoretical work has been done on QSL states in 
systems with strong spin orbit coupHng, i.e. with strongly bro- 
ken SU(2) symmetry. As such, the structure of possible QSL 
states in the present model requires particular investigation. 
Guided by the pseudo-singlet structure, and the hidden SU(2) 
symmetry of the model, we construct candidate QSL states for 
the full Hamiltonian, Hex-i, by a slave-particle technique. 



A. Spin waves 

In the previous two sections, we have discussed the state 
phase diagram based on mean-field theory. Here, we perform 
a linear spin wave analysis, which perturbatively describes the 
effect of quantum fluctuations on the various phases obtained 
so far, and also predicts the structure of collective modes, 
which might, e.g., be observed in inelastic neutron scattering. 
Finally, because we have not explored the full space of mean- 
field states, the calculation also provides an important check 
that the phases we have found are at least metastable. 

The conventional Holstein-Primakoff (HP) transformation 
for spin-S operators cannot be directly applied for the three 
variational ground states because none of the three states is 
a fully polarized state for any projection of the spin angu- 
lar momentum operator j. This is especially severe for the 
AFM state, for which the spin expectation value simply van- 
ishes. Instead, we formulate an "SU(4) spin wave theory", by 
rewriting the Hamiltonian, ( [27] ), in a bilinear form in terms 
of the 15 generators of the SU(4) group. To do so, we in- 
troduce, for any local basis for the single-site Hilbert space 
{|n)}, n = 1, 2, 3, 4, the complete set of operators^^ 

SZ = \m){n\. (85) 

These SU(4) generators obey the algebra [SJ^^Sl] = 
^nkSly^ — SmiS]^. We can then use the HP transformation for 
the generators of SU(4). In this transformation, one selects a 
particular state in the four-dimensional basis to be the vacuum 
and introduces three bosons associated with excitations to the 
three other states. For the AFM phase, the classical ground 
state is given by Eq. ( |40| ). Hence we take, on the A sublattice 
(zi integer), the basis 
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The Hamiltonian H in Eq. ( |27l > in this basis has a quadratic 
form. 



(88) 



where the coefficients Ckimnihj) (which are straightforward 
to obtain, so we do not give them explicitly) depend linearly 
on J, J' and V. 

To introduce the HP transformation on sublattice A (B), 
we choose |1) as the vacuum which is annihilated by three 
"magnon" annihilation operators (6^), n = 2,3,4. One 
can think these three bosons as descending from mixed spin 

and orbitals fluctuations of the Hamiltonian before the V3 

2 

projection. The HP transformation is defined as^^, for i in 
the A sublattice. 



Sl{i)=M- 



^a],(z)a^(z). 



(89) 



^n(0=4(^)J^-E^K^K«, (n^l) (90) 

V 19^1 

Sl{i)=ai{i)a,{z). (/,n7^1),(91) 

while for i in the B sublattice, the same formula holds with 
an{i) replaced by bn{i). In the above equations M is defined 
as the number of columns in the Young tableaux for the repre- 
sentation of SU(4). In our case (fundamental representation), 
we must set M = 1. In the generalization to arbitrary M, the 
classical limit where the classical ground state becomes ex- 
act is M 00. However, we apply this HP transformation 
directly for M = 1. Inserting this into ( [88] ), we expand it to 
obtain a quadratic form in the bosonic operators. The constant 
term in the expansion gives the classical ground state energy. 
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independent of the phase 0. The quadratic terms lead to quan- 
tum corrections. Defining the Fourier transform of the bosonic 
operators, the spin- wave Hamiltonian can be organized in the 
form i^k with 



A 



A 



(93) 



where A = («2k, «3k, <^4k, ^2k, ^3k, ^4k) is the vector of 
magnon annihilation operators and Fk and Gk are 6 x 6 ma- 
trices. This spin- wave Hamiltonian can be diagonalized by 
standard methods. 

For the AFM ground state, Eq. ( |4Q|), we obtain a gapless 
"magnon" mode, as depicted in Fig. |14[ This gapless mode 
is associated with the continuous accidental degeneracy, and 
indeed occurs for arbitrary (j). By contrast, in the FMllO and 
FMIOO phases, one observes a gap for all the spin wave modes 
(see Fig. [15] and Fig. [16]). The gap in the FMllO phase in- 
creases with J', as expected since this corresponds to increas- 
ingly violated SU(2) symmetry. In all cases, the modes are 
all well-defined with positive real frequencies, indicating the 
stability of the phases in the classical sense: i.e. that we have 
properly found local energy minima of the mean field theory. 

Finally, having obtained the spin wave modes, we can eval- 
uate the quantum corrections. It is most interesting to consider 
the reduction of the order parameter by quantum fluctuations. 
We can define this by considering the probability to find a 
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given spin in its mean-field ground state. This is nothing but 
the vacuum state of the HP bosons. Hence this probabiHty is 
given, for a site on the A sublattice, by 



(94) 



This quantity is directly analogous to the staggered mag- 
netization in the usual HP treatment of a quantum antifer- 
romagnet. We therefore denote AM = 1 — Pgs{i) = 
'}2n^ii^ni^)^n{^))- ^^^^ obtained, at T = 0, by integrat- 
ing the zero point contribution to the boson number from each 
spin wave mode. The quadratic spin- wave Hamiltonian is di- 
agonalized by a Bogoliubov transformation Qk, 



(95) 



in which, Ck = (cik, C2k, csk, C4k, C5k, cek), and Qk is a 12 x 
12 matrix. From this we obtain the quantum correction 
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Numerically, we find that this quantum correction is maxi- 
mal for J' = V = 0, and is given by AM ^ 1.7 at this point. 
This is much larger than 1, implying that the fluctuations at 
this point are large and that the mean field theory is at least 
quantitatively invalid. For increasing J' and V the correction 
becomes significantly smaller, and mean field theory may be 
reliable. In the vicinity of the J' = F = 0, one may expect 
a very different ground state, incorporating strong quantum 
fluctuations. We explore some possible non-magnetic ground 
states in the remainder of this section. 




FIG. 14. Spin wave spectrum for the AFM phase at J' = 0.1 and 
y = 0.2 along [001] momentum direction. And J — 1. There is one 
low-lying gapless mode. The fee lattice constant is set to be a = 1. 
And the phase = for the ground state in Eq. (Hol. 



B. Pseudo-singlets in different planes 

We start our analysis of non-magnetic states by considering 
two sites in the XY plane, which interact with the Hamiltonian 




FIG. 15. Spin wave spectrum for the FMllO phase at = 0.3 and 
V — 0.2 along [001] momentum direction. The lowest excitation 
mode has an energy gap A = 0.241 at /c = 0. In the graph, J — 1. 
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FIG. 16. Spin wave spectrum for the FMIOO phase at = 0.1 and 
V — OA along [001] momentum direction. The lowest excitation 
mode has an energy gap A = 0.0224 at /c = 0. In the graph, J = 1. 



^eS^i i)- Remarkably, the ground state has a form identical 
to an SU(2) spin singlet, if written in terms of pseudospin-1/2 
states = ±1/2: 



2)3 I 2)^12)^) 



(97) 



One may understand this result by writing down the projected 
spin and occupation number operators in xy orbitals, in the 
basis of eigenstates (see Eq.[23]): 
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(98) 



(99) 



0/ 



Here cr is the vector of Pauli matrices and I2 is the 2x2 iden- 
tity matrix. One may consider Sxy as an effective spin- 1/2 
operator in the subspace of = ±1/2 states, which naturally 
explains the SU(2) singlet in Eq. ([97]). 

For the XZ and YZ planes, one simply needs to apply a cu- 
bic permutation to the results obtained for XY planes, or more 
formally, apply a unitary transformation that rotates about the 
[111] axis by ±27r/3. 
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FIG. 17. (color online) Columnar Valence Bond Solid (VBS) state 
within an XY plane. The dashed square indicates the face of a con- 
ventional cubic unit cell, while the solid lines connect the FCC near- 
est neighbors within the plane, which form a 45° rotated square lat- 
tice. 



with 



U = exp 



' 3 



(102) 



The upper indices jj.' = p{ii) and ji" = p~^{lj) in 
Eqs. ( |1Q0| ) and ( |101| ) denote cyclic and anti-cyclic permuta- 
tions of x^y^z, respectively [i.e. p : {x^y^z) {y^z^x), 
with inverse : (x, y, z) (z, x, y)]. The two-site ground 
states in the XZ and YZ planes are the pseudo-singlets in the 
subspace of = ±1/2 states and = ±1/2 states, respec- 
tively. 



C. Valence bond solid state 

It is natural to consider a product state of such pseudo- 
singlet "valence bonds" (also called "dimers") as a candidate 
(prototypical variational) non-magnetic ground state. To do 
so, we must divide the spins into two neighboring sublat- 
tices, which will be paired. This by necessity breaks lattice 
symmetries. Such a state is called a Valence Bond Solid, or 
VBS, state. At the level of valence bond product states, many 
possible arrangements of the dimers are degenerate. This de- 
generacy is artificial and will be broken if the wavefunctions 
are improved. We will not investigate this in any detail, and 
just consider the simplest VBS state in which the dimers form 
a "columnar" arrangement within a single (001) plane. See 
Fig. 
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The variational energy of such a state (actually any state 
with a planar arrangement of dimers has the same energy) is 
readily evaluated. We obtain the energy per site EyBs/J^ = 
{VBS\Hex-i\VBS)/M = -5/12J = -0.42J. This is 
slightly higher than the mean-field ground state energy of the 
AFM state, E^p^jM = - J/2. However, the large quantum 
fluctuations are expected to destabilize the latter state, and per- 
haps might stabilize the VBS one. So such a VBS state seems 



competitive, and may be considered as a possibility for future 
exploration. 



D. QSLs and Fermionic mean field theory 

The most general approach that has been applied to de- 
scribe QSL states is the slave particle method, in which aux- 
iliary fermions are introduced, and the ground state for the 
spin system is described by some projection of a nominally 
simple fermionic state into the physical spin Hilbert space. 
This results, in the usual SU(2)-invariant case, in wavefunc- 
tions which are composed of superpositions of products of 
SU(2) singlets. Here, the appearance of two-site pseudo- 
singlet ground states points to the possibility of applying a 
similar fermionic mean field theory. In this section, we imple- 
ment this technique for the full antiferromagnetic exchange 
interaction. 

We first introduce the auxiliary fermionic creation opera- 
tors, whose quanta we call "spinons'l^ 



= /iLk^cuum), a = 1, ...4, 



(103) 



where for convenience we have relabeled the states = 
|,^,— |bya = 1,2, 3, 4, respectively. The physical 
Hilbert space is constructed from states with one fermion at 
each site, which imposes the constraint 
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In this notation, the spin and number operators become 



S F 

F^rf- F. 



(105) 
(106) 



where on the right-hand side it is to be understood that the 
matrices in Eqs. ([98]) and ([99]) act on the vector of spinon op- 
erators 



(/il, /i2, /i3, fi4.Y 



(107) 



Similar expressions can readily be written for operators in XZ 
and YZ planes. Thus the Hamiltonian in terms of spinons 
reads 



2J 
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q;,/3=2,3 



(108) 



in which, ctq,^ • a^^fd' = 25^(3' S a' (3 — ^ajS^a'jS' has been used. 
Similar spinon Hamiltonians can also be written down for XZ 
and YZ planes. When we write down the full antiferromag- 
netic exchange Hamiltonian and sum over XY, YZ and XZ 
planes, we find that the single- site terms, which are quadratic 
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in spinon operators, sum up to a constant once we impose the 
single occupancy constraint. We are then left with the terms 
that are quartic in spinon operators. 

We now follow the standard procedure of slave particle 
mean field theory to decouple the quartic terms in the spinon 
Hamiltonian and write down a mean field ansatz. We start 
with the exchange Hamiltonian in the XY plane, Hf^.i . 

We require a mean field ansatz for the fermionic bond ex- 
pectation values, {fiafj(3)' Noting the structure of the two 
site pseudo-singlet in this plane, we choose an ansatz which 
reproduces a quantum ground state of this type. Specifically, 



(fiafjs) ~ Xiji^xy)(: 



with 



















0/ 



(ij) G XY (109) 



(110) 



Note that, by construction, this expectation value is invariant 
under the hidden SU(2) symmetry. The Xij on the XZ and YZ 
planes are determined by symmetry 



and 



Xij;al3 = ifLfjfi) = Xij{1xz)l3a, {ij) & XZ 
Xii;a/3 = ifLfils) = Xii{Ivz)/3a, {ij) G YZ 



(111) 



(112) 
(113) 



with the unitary transformation introduced in Eq. ( [1Q2| ). 
We then arrive at the mean field Hamiltonian 



XY 
MF 



{ij)exY 



(114) 



with J = 2 J/9. Here are the Lagrange multipliers related 
to the single-occupancy constraint in Eq. (|1Q4| ). ^mf 
can be readily written down using Eq.|lll and Eq. (|112|). 



1. uniform spin liquid 

As discussed in Sec. the antiferromagnetic ex- 
change Hamiltonian has a "hidden" global SU(2) symmetry, 
[G^,?iCex-i] =0. It is easy to find that the full mean field 
Hamiltonian we have here respects this "hidden" global SU(2) 
symmetry. We seek a quantum spin liquid ground state which 
does not break any symmetries of the original Hamiltonian. 
Translational invariance imposes A^ = A = const. First we 
consider the ansatz for a uniform spin liquid 



■ Xji 



for i, j nearest neighbors on the fee lattice. This naturally re- 
spects point group and time reversal symmetries. The Hamil- 
tonian in Eq. ( [114| ) is then diagonalized by Fourier transform 



/„(k) = ^e-'^-^Via. 



(116) 



We find 

nMF= E E^^(k)[/t^(k)/,+(k) + /t_(k)/,_(k)] 



A=l,2 k 



+12A/'JX^ (117) 
where A = 1 , 2 label doubly degenerate bands with dispersion 



ei,2(k) = Jxei,2(k) = -2Jx 



CxCy 



- CyCz 



■ CzCx 



^x^y ~^ ^y^z ~^ ^z^x ~ ^xCyC z{Cx-\-Cy-\-C z) 

(118) 

Here Cx = cos{kx/2), Cy = cos{ky/2) and Cz = 
cos{kz/2). The double degeneracy of the two bands is due 
to Kramer's degeneracy, since ji is a spin-3/2 operator and 
the Hamiltonian has time reversal symmetry. 

The ground state wave function at the mean field level is 
described by a Fermi sea of spinons 

|*MF>= n n/l+(k)/l-(k)|vacuum), (119) 

A=l,2 k 

for all k below the Fermi surface. The mean field ground state 
energy per site is 



J\fJ 



/(^%-e~A(k)]6,(k) + l2x^ 



A=l,2 



(120) 

where the integral is over the first Brillouin zone of the fee 
lattice and the dimensionless chemical potential /i is fixed by 
the quarter filling condition 



2E /7S^^['^-^^^(k)] = i- 

A=l,2'^ ^ ^ 



(121) 



The wave function in Eq. \\\9\ must be Gutz wilier-projected 
into the physical Hilbert space with one spinon per site 



1^) =K.=i|^mf). 



(122) 



Here we simply evaluate the ground state energy at the 
mean field level. Minimizing Eq. ( |12Q| ) with respect to the 
parameter x, we find 



X 



\ 1 o 



A=l,2 ' 



d^h 

^^[/i-6A(k)]6A(k), 



and 



E 



MF 



= -lJ{xl'- 



(123) 



(124) 



The spinon density is at quarter filling for fi ^ —1.58 J. The 
mean field energy for the uniform spin liquid state is then 

n(0) 



(115) £;]^^^/A/' ^ -0.041J. 
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2. 7T-flux spin liquid 



We now consider the ansatz for the tt flux spin liquid state 



iUustrated in Fig. 18 In order to preserve time reversal sym- 
metry, the phase of the Xij at each bond can only assume the 
values or TT. We divide the fee lattice into four cubic lattices 



'J' 



r/ = (0,0,0)+R 
rf = (1/2, 1/2,0) + R 
= (0,1/2, 1/2) + R, 
rf = (1/2, 0,1/2) + R, 



■J' 



where Hj is a unit vector in the cubic lattice with lattice pa- 
rameter a = 1. We denote by Fja the vector of spinon an- 
nihilation operators at site j of sublattice A, and similarly for 
the other sublattices. We assign Xij = +X ^ the bonds con- 
necting sites in sublattice A to all its nearest neighbors in sub- 
lattices B^CjD, and Xij = — X ^ the bonds connecting two 
sites that belong to sublattices B^C or D. In other words, 
this ansatz corresponds to assigning — x to the three bonds in 
the BCD base of each tetrahedron in the fee lattice and +x 
to the three bonds connecting the BCD base to the A vertex. 
As a result, there is tt flux through every triangle and zero 
flux through every square in the fee lattice. While this ansatz 
is clearly invariant under point group transformations about 
A sites, it is also invariant under lattice translations, despite 
the fact that this permutes the 4 sublattices. This is because 
the corresponding changes in Xij can be removed by a gauge 
transformation. For instance, the gauge transformation 



Fjc ^ F, 



(125) 



exchanges the signs of Xij between sublattices A and B. It 
follows that the tt flux ansatz is invariant under point group 
symmetries about any site of the fee lattice and therefore re- 
spects all symmetries of the original Hamiltonian. 

Minimizing the energy for the mean field Hamiltonian with 
four sublattices, we find 8 doubly degenerate bands for the tt 
flux state. Quarter filling is reached for dimensionless chemi- 
cal potential ji ^ —1.68 J. The mean field energy is 



Emf 
J\fJ 



8 

A=l 



(2^ 



^[/i-eA(k)]eA(k) + 12x', (126) 



where the integral is over the reduced Brillouin zone of the cu- 
bic lattice. We find E^^\/J\f = -0.053 J. This is lower than 
the energy for the uniform state. We may also compare this to 
the mean field energy of a VBS state (it is not so meaningful 
to compare these slave particle mean field energies directly 
to the variational energies quoted earlier for the Weiss mean 
field and VBS states). For a mean-field VBS state, we take Xij 
non-zero only on a set of non-overlapping dimers. In this case, 
we obtain E^^g/M = -1/18J ^ -0.055 J. This is slightly 
lower than the QSL states, but we expect that the energy of the 
spin liquid states will be lowered by the Gutzwiller projection, 
since the latter is known to enhance spin-spin correlations.!^ 




FIG. 18. (Color online) Conventional unit cell of the fee lattiee di- 
vided into four sublattices. The ansatz for the tt flux state corre- 
sponds to assigning hopping amplitude +x to the bonds represented 
by solid lines and — x to the bonds represented by the dashed lines. 



We note that both spin liquid states have Fermi surfaces 
which are not nested, and have no obvious instabilities. The 
states are also stable against bond anisotropy which enhances 
the hopping in a given plane. Perturbations to NN AFM 
exchange, such as next-nearest-neighbor interactions, will in 
general require more general ansatze for the bond matrix Xij • 
Nonetheless, as long as the perturbations are in some sense 
small, the Xij assumed in Eq. ( |109| ), in which hopping in a 
given plane occurs predominantly for two out of four spinon 
species, should be a good starting point for approximations. 



VI. DISCUSSION 

In this paper, we have introduced and analyzed a model 
to describe localized electrons in a 4d^ or 5d^ configuration 
on an fee lattice, in which strong spin-orbit coupling and 
the t2g orbital degeneracy combine to produce an effective 
j = 3/2 description. The model contains three interactions 
- nearest neighbor antiferromagnetic and ferromagnetic ex- 
change, and electric quadrupolar interactions - and in addi- 
tion may include the effect of structural anisotropy. We obtain 
the (Weiss) mean-field phase diagram, which includes 3 main 
phases, which all have a two- sublattice Q = 27r(001) struc- 
ture. In all the phases, large multipolar order parameters in 
addition to the usual magnetic dipolar order are present. Most 
remarkably, we find a broad regime of time-reversal invariant 
but quadrupolar ordered phase at intermediate temperatures. 
A spin-wave analysis indicates that quantum fluctuations are 
strong when nearest-neighbor antiferromagnetic exchange is 
dominant, and in this case, we suggest possible quantum spin 
liquid and valence bond solid phases. 
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A. Experimental ramifications 

The theory developed here can be applied and tested in a 
multitude of ways. Here we discuss a few of the main exper- 
imental properties which might be measured. First, there is 
the spatial symmetry breaking of the ordered phases. All the 
ordered states, at least in the cubic case, break lattice symme- 
tries, and in particular double the unit cell. This can be tested 
in experiments such as neutron and x-ray scattering. 

We discuss in some further detail the most intriguing case of 
the quadrupolar phase, which is non-magnetic. It constitutes a 
type of real (time-reversal invariant) orbital ordering. It leads 
to a spontaneous splitting of the local quadruplet, breaking it 
down to an elemental Kramer's doublet. As such, this is not 
entirely distinct from a Jahn-Teller effect, in which ionic mo- 
tions would lead to such splitting. In particular, even though 
in our model atomic displacements are not involved in an es- 
sential way, they would be expected at least to respond to the 
orbital order. In principle, this could be measured by scatter- 
ing (x-rays or neutrons) which accurately measure the crystal 
structure and symmetry. From the order parameter description 
of the quadrupolar phase, we can obtain the corresponding 
space group and crystal structure parameters to be sought in 
such a measurement. In particular, we find that the quadrupo- 
lar ordered phase corresponds to the tetragonal space group 
P42/mnm (number 136). In this space group, apart from the 
doubling of the unit cell, one finds that all the A sites, B sites, 
and B' sites remain equivalent. However, the oxygens are no 
longer equivalent, but split into three classes, occupying the 
4e, 4f, and 4g Wyckoff positions. Each of these positions 
has one degree of freedom which is not fixed by symmetry. 
Physically, the oxygens remain constrained to the simple cu- 
bic axes of the perovskite reference unit cell, but may move 
by different amounts along each of the three axes. This is two 
more degrees of freedom than is found in the cubic Fm3m 
(number 225) space group, in which the oxygens maintain an 
ideal octahedron equidistant from each B (or B') site. While 
symmetry requires these structural changes, we do not have at 
present an estimate for their magnitude, which could be weak 
if coupling to the lattice is not strong. 

Another quantity we have already discussed in Sec 
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is the magnetic susceptibility, which shows signatures of the 
quadrupolar and ordering transitions. One complication is that 
the susceptibility is in many cases highly anisotropic, and one 
may not be sure what component(s) are being measured in 
practice. Specifically, one may expect, if the system is ideal 
and fully in equilibrium, that the broken symmetry order pa- 
rameters can be reoriented by the magnetic field, in such a 
way that they minimize the free energy. This will typically fa- 
vor orientations which maximize the magnetic susceptibility. 
For instance, in the quadrupolar phase for cubic symmetry, 
this is an orientation with Q perpendicular to the field. How- 
ever, such reorientation involves motion of domain walls and 
very large numbers of spins, and can easily be prevented by 
pinning or at least be incomplete in practice. Thus some di- 
versity of behavior may be expected in experiment, as well as 
possibly hysteretic behavior even in the non-magnetic state. 
When the crystal is non-cubic, one may explore the influ- 



ence of single-ion anisotropy on the magnetic susceptibility. 
A naiVe application of Eq. ([84]) would immediately imply that 
single-ion terms do not contribute to the Curie- Weiss tem- 
perature as measured in the powder susceptibility. However, 
we caution that these equations hold only in the true high- 
temperature regime, in which T ^ \D\. If T is smaller than or 
comparable to higher order terms in the high temperature 
expansion are non-trivial, and a non-vanishing fitted Curie- 
Weiss temperature may result from D alone. Let us consider 
the powder susceptibility x for independent ions (i.e. neglect- 
ing exchange). One has 



X 



Xzz 

~3~ 



2X. 
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4T 



tanhL)/T 
2D ' 



(127) 



in units of g^jJ^^. We then suppose a linear fit to 1/x ver- 
sus T is made over a narrow region in the neighborhood of 
the temperature T^t, and extrapolated to find the Curie-Weiss 
temperature as the intercept of the horizontal axis. The result 
is 



©cw(^fit) 



2Tfit (Tfit sinh(|£)-2I)) 
L>(3cosh(|£) + 7) 



(128) 



Note that the fitted Curie- Weiss temperature is always nega- 
tive, and is independent of the sign of D. It reaches a max- 
imum in magnitude (at fixed D) of 9cw ^ — 0.18|I)| when 
Tfit ~ 0.88 1 T) I , and only approaches zero very slowly when 
Tfit ^ \D\: 6cw ~ Co^v^^^s^IY' ^ fixed fitting 



temperature, the maximum achievable Curie- Weiss tempera- 
ture is Gcw ^ -0.26Tfit, when \D\ = 1.84Tfit. 



B. Materials survey 

We now turn to a discussion of specific materials which 
have been studied experimentally. 



1. Ba2YMoOQ 

We begin with the material Ba2YMo06, which has been 
suggested experimentally to be an exotic "valence bond glass" 
or to have a "collective spin singlet" ground state. The ex- 
pected separation between the j = 3/2 and j = 1/2 states 
in this material is over 2000K, so that the effective j = 3/2 
description used here should be excellent. Two recent exper- 
imental papers^^ observed an unusual behavior of the mag- 
netic susceptibility, with two Curie regimes, such that 1/x is 
linear both above lOOK and below 50K. Moreover, the mag- 
netic specific heat shows a peak around 5 OK, with Ref. |36] 
estimating the total magnetic entropy approximately equal to 
Rln4, as expected for j = 3/2. Both these results suggest 
the existence of some single-ion anisotropy, which would ex- 
plain the existence of two Curie regimes because it splits the 
4-fold degeneracy of the j = 3/2 states but leaves a 2-fold 
Kramer's doublet at temperatures below \D\, which still gives 
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a Curie signal. However, the cubic symmetry observed experi- 
mentally seems to rule out such an explanation. Moreover, the 
form of the powder susceptibility in Refs. 36|and[37]is qualita- 
tively different from that expected for either fixed easy-plane 
or easy-axis anisotropy. 

These difficulties are resolved if one considers the possibil- 
ity of spontaneous anisotropy, which indeed is the primary 
characteristic of the quadrupolar ordered state. For exam- 
ple, the mean-field susceptibility for the cubic model with 
J' = 0.2 J, V = 0.3 J is plotted in Fig. fTsl At tempera- 



tures above the FMllO phase, one indeed observes two Curie 
regimes in the susceptibility, with a larger Curie constant at 



low temperature, as seen in the experiments. The kink in x co- 
incides with the quadrupolar ordering transition, and there is a 
peak in the specific heat at this temperature, also as observed 
in experiment. The theoretical specific heat has a second peak 
at lower temperatures, associated with magnetic ordering and 
exhaustion of the unsplit Kramer's doublet. We suggest that 
this peak is below the lowest temperatures measured, or per- 
haps is avoided due to disorder, and the spins falling out of 
equilibrium at low temperature. The fact that the Curie-Weiss 
temperature extracted below 50K is only -2.3K corroborates 
the notion that any magnetic ordering may be too low to ob- 
serve or be obscured by the effects of disorder. 



Compound 


config. 


crystal structure 


Gcw 


Meff(Ms) 


magnetic transition 


frustration parameter / 


Ref 


BasYMoOe 


Mo^+(4(i') 


cubic 


-91K 


1.34 


PM down to 2K 


/>45 


I38J 


BasYMoOe 


Mo^+(4d') 


cubic 


-160K 


1.40 


PM down to 2K 


/ ^ 80 


1361 


Ba2YMo06 


Mo^+(4(i') 


cubic 


-219K 


1.72 


PM down to 2K 


/> 100 


137] 


La2LiMo06 


Mo^+ (4d') 


monoclinic 


-45K 


1.42 


PM to 2K 


/ ^ 20 


137] 


SrsMgReOe 


Re^+(5(i^) 


tetragonal 


-426K 


1.72 


spin glass, To 50K 


/>8 


|39j 


Sr2CaRe06 


Re^+(5d') 


monoclinic 


-443K 


1.659 


spin glass, Tg ^ 14K 


/>30 


140] 


Ba2CaRe06 


Re^+(5(i') 


cubic to tetragonal (at T 120K) 


-38.8K 


0.744 


AFMTn = 15.4K 


/~2 


HH 


Ba2LiOs06 


Os^+(5d^) 


cubic 


-40.48K 


0.733 


AFM Tn ^ 8K 


/>5 


|42| 


Ba2NaOs06 




cubic 


-32.45K 


0.677 


FM Tn - 8K 


/>4 


142] 


Ba2NaOs06 


Os^+(5(i^) 


cubic 


lOK 


- 0.6 


FM Tn 6.8K 


/>4 


114] 



TABLE II. A list of ordered double perovskites. Note the discrepancy in Curie-Weiss temperature and /ieff may originate from the experimental 
fitting of data at different temperature range. 



2. La2LiMoOQ 



3. Sr2CaReOQ and Sr2MgReOQ 



La2LiMo06 is monoclinic, the deviation from cubic sym- 
metry arising primarily from rotations of the octahedra. The 
local coordination of the Mo sites is nearly perfectly octahe- 
dral with a weak tetragonal compression. The nature of crystal 
field effects, if significant, is unclear at present. Magnetically, 
the susceptibility shows, like Ba2YMo06, two apparent Curie 
regimes, separated by a kink at approximately 150K. How- 
ever, opposite to that material, La2LiMo06 shows a smaller 
effective moment at low temperature compared to high tem- 
perature. In addition, the high temperature Curie-Weiss tem- 
perature is 9cw ^ — 45i^r, significantly smaller than the kink 
temperature. Irreversibility distinguishing the behavior of the 
ZFC/FC susceptibility appears below 25 K. 

The appearance of two Curie regimes again suggests either 
fixed or spontaneous magnetic anisotropy setting in around 
150K. However, the reduction of the effective moment below 
the kink in is puzzling. We did not find this behavior in 
the powder susceptibility within our model, with or without 
anisotropy modeled by D. As remarked above, however, the 
actual nature of the crystal field anisotropy in La2LiMo06 is 
unclear. If it is significant and different in form from the D 
term, this might explain the behavior. Single crystal studies 
would be helpful in elucidating the situation. 



Sr2CaRe06 and Sr2MgRe06 have distorted perovskite 
structures, with monoclinic and tetragonal symmetry, 
respectively.^^ Experimentally, the materials are notable 
for their very high antiferromagnetic Curie- Weiss tempera- 
ture, — 6cw ^ 400K. Susceptibility and specific heat mea- 
surements show anomalies suggestive of freezing and/or 
short-range ordering at 14K and 50K, for Sr2CaRe06 and 
Sr2MgRe06, respectively. Two possible interpretations of this 
behavior are: (1) the Curie- Weiss temperature is dominated by 
strong exchange, but fluctuations largely suppress ordering, or 
(2) the Curie-Weiss temperature is due largely to single-ion 
effects, and the true exchange scale is comparable to the ob- 
served anomalies in x and Cy . 

In the former scenario, the key question is why these two 
materials show so much larger exchange than do the other 
compounds in this family. Erom the point of view of this 
work, attributing the Curie- Weiss temperature to exchange 
alone would imply J is actually comparable to the SOC, so 
that the projection to j = 3/2 may even be suspect. The 
Curie-Weiss temperatures are sufficiently large that one may 
suspect that the bd electrons are not so well localized, and the 
system is close to a Mott transition. It would be interesting to 
measure their optical properties to address this possibility. 

The latter explanation seems possible, as both materi- 
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als show significant deviations from the cubic structure: 
Sr2CaRe06 is monocHnic, while Sr2MgRe06 is tetragonal. 
The actual distortions of the octahedra are rather small in both 
cases, the Re-0 distance varying by only about 0.02A at room 
temperature. However, there are significant rotations and tilts 
of the octahedra, and crystal field splittings of the j = 3/2 
quadruplet are certainly allowed. Examination of the Re-0 
bond lengths suggests easy-axis anisotropy. From Eq. ( |128| ), 
we see that in principle a negative Curie-Weiss temperature 
could be attributed to D. However, from the present model we 
cannot obtain such a large value, which in these two materials 
is comparable or larger than the fitting temperature. Neverthe- 
less, we may imagine that some combination of exchange and 
single-ion anisotropy may conspire to produce the observed 
behavior. 

If we assume a large easy-axis anisotropy, we would then 



expect, based on the the analysis in Sec. Ill A 1 to have an 
AFM ground state. The anomalies might be related to this 
ordering. Experimentally, spin freezing and irreversibility is 
observed, but without clear signs of long-range ordering. The 
experimentalists caution that, due to the small magnetic mo- 
ment of the Re^+ ions, a sm all ord ered component could not 
be ruled out in either material .^^^l^ Indeed, in the AFM state, a 
very small moment is expected, due to the primacy of octupo- 
lar order. 

While this is promising, we note that it is likely that sev- 
eral effects not in our model play a role. First, the struc- 
ture of the materials is not a simple compression of the cu- 
bic structure, and so the crystal fields might have a signifi- 
cantly different form from the simple D term. This is espe- 
cially true in Sr2CaRe06, which has the more distorted mon- 
oclinic structure. Second, in Sr2MgRe06, the Re-O-Mg bond 
angles are very different in the XY plane (160°) and normal 
to it (180°),^^ so substantial spatial anisotropy in the exchange 
couplings may be present. This is not included in our model. 
Also in Sr2MgRe06, the ZFC and EC susceptibility actually 
diverge already around 300K, which suggests a high degree 
of disorder in this material, which might be responsible for 
converting the AFM to a glassy state. 

One indication supporting an antiferromagnetic ground 
state is the observation, in Sr2CaRe06, of a magnetic 
contribution to the specific heat, in contrast to the usual lin- 
ear one characteristic of a spin glass. The behavior 
would naturally be expected from spin waves in the AFM 
state, which as we have noted displays gapless spin waves, at 
least in the semi-classical approximation. Such behavior 
might even persist if the AFM order had a finite correlation 
length, due to Halperin-Saslow modeJ^, as recently postu- 
lated in NiGa2S4E2. A T-linear specific heat was observed 
in Sr2MgRe06, albeit with a small coefficient.^^ As we have 
already remarked, however, this material is likely to be more 
disordered, consistent with the more conventional spin glass- 
like specific heat. 

While these considerations seem reasonable, they are 
hardly definitive. Further studies, particularly on single crys- 
tals, would be most helpful in clarifying the physics of these 
materials. 



4. Ba2CaReOQ 

In Ba2CaRe06, there is a structural transition from a high 
temperature cubic phase to a low temperature tetragonal one, 
with a doubled unit cell, at T = 120K'^^ The experimen- 
talists have fitted the low temperature structure to the I4/m 
space group. From this fit, they found an elongation along 
the c (or z) axis, but a slight compression of the ReOe oc- 
tahedra. One may consider two possibilities. Either this is 
indeed the correct symmetry, in which case it must have struc- 
tural origin not related to the 5d electrons, or this transition in 
fact coincides with the quadrupolar ordering described here, 
which also gives a tetragonal state with the same unit cell. 
The P42/mnm space group was not considered in Ref 41 

A small, negative Curie- Weiss temperature Bcw = — 39K 
was measured by fitting the susceptibility in the high tempera- 
ture cubic phase , indicating that J' should be not too large. A 
predominantly antiferromagnetic ordering transition was ob- 
served at T =15K, which is consistent with our expectations 
in the small J' and V regime (recall that an antiferromagnetic 
state is expected both with and without single ion anisotropy 
in this parameter range). It would be interesting to compare 
the predicted magnetic structure in the AFM or AFM' phase 
with experiment, by carrying out neutron scattering and NMR 
measurements. 



5. Ba2NaOsOQ 

Ba2NaOs06!^is one of only two examples in this class in 
which single crystal experiments have been performed, to our 
knowledge. A transition at T =6.8K was found to a ferromag- 
netic state with easy axis along a [110] direction. Within the 
experimental resolution, the material was found to remain cu- 
bic down to the lowest measured temperature. A fit to Curie- 
Weiss behavior found a negative Curie- Weiss temperature, 
— 13K< ©cw < — lOK, depending upon field orientation. We 
note that this is compatible with Eq. ([84]), and suggests that the 
material is in the regime of larger V and smaller J^ In this 
region, we expect a high temperature quadrupolar transition, 
well above the ferromagnetic state. While no such transition 
is observed in the experiment, the magnetic specific heat of 
such a transition may be masked by the lattice contribution at 
higher temperature, and the signature in susceptibility may be 
subtle. 

Several other indications are in favor of this scenario. First, 
the susceptibility continues to display pronounced anisotropy, 
favoring the [110] direction, up to at least 200K, well above 
the ferromagnetic transition. This would indeed be expected 
in the quadrupolar ordered state. Second, the magnetic spe- 
cific heat divided by temperature, integrated over the peak 
up to about 15K gives a magnetic entropy of approximately 

In 2, which is only half the expected entropy for the j = 3/2 
quadruplet. This entropy may be released over a signifi- 
cantly higher temperature range, up to the quadrupolar or- 
dering transition. Third, the observation itself of ferromag- 
netism with a [110] easy axis is a marked success of our the- 
ory. This type of anisotropy is not natural from the standpoint 
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of Landau theory, within which cubic anisotropy manifests it- 
self at leading order as a term in the free energy of the form 
v[(M^)'^ + (M^)'^ + (M^)'^] (M is the magnetization), which, 
depending upon the sign of generates a [111] or [100] easy 
axis. A continuous Landau transition to the FMllO state is 
instead made possible by the fact that the quadrupolar order 
already breaks the cubic symmetry in the paramagnetic state. 

It would be interesting to further probe the system to estab- 
lish in more detail the correspondence (or lack thereof) with 
our theoretical predictions. The predicted tetragonal distor- 
tion of the cubic structure would be a natural quantity to seek 
in experiment. This also manifests locally in the magnetism, 
since although the net magnetization is aligned with the [110] 
axis, the local spin expectation values are not. This might be 
measurable for instance by a local probe such as NMR. 



6. Ba2LiOsOQ 



Ba2LiOs06 has also been grown in single crystal form. 
The structure was determined to be cubic by x-ray diffrac- 
tion at room temperature."^^ Aside from this diffraction data, 
only bulk magnetic susceptibility results are available. One 
observes a negative Curie- Weiss temperature 6cw ^ — 40i^, 
and an apparent antiferromagnetic transition at Tn ^ 8i^. 
This appears largely consistent with the expected behavior in 
region I of Fig. [6] 



C. Comparison and future work 

We are aware of only one other theoretical work studying 
this class of materials. Lee and Pickett^^ performed electronic 
structure calculations for Ba2NaOs06 and Ba2LiOs06, em- 
phasizing the role of SOC. We completely agree with the con- 
clusion that SOC plays a crucial role in the magnetism. How- 
ever, the magnetic structure and phase transitions were not 
addressed. 

Our study is much more comprehensive, and gives a great 
deal of guidance both for future theory and experiment. Many 
experimental suggestions have already been made. In partic- 
ular verification of the quadrupolar ordering transition would 
be especially exciting. On the theoretical side, the problem of 
the effects of quantum fluctuations in the small J' and V limit 
remains rather open. It would be remarkable if a spin liquid 
or valence bond solid state could be established for this highly 
non-SU(2) symmetric and nominally "large" spin j = 3/2 
model. To do so would require some hard theoretical work 
applying more quantitative numerical methods to our model 
Hamiltonian. A natural extension of this work would be to 
consider the "higher spin" analogs of these materials, with 
4d^ or 5d^ electronic states. As there is still partial occupation 
of the t2g orbitals in this case, we expect SOC again to play 
a dominant role, and interesting multipolar physics is likely 
present. 
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